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r I m t^iL Z ^Tr. ^ a PP ,,ca, ? r us f! ada P tive nullin 8 »n<l focusing with non-invasive electric field probes 
(1 2) to control the electnc field intensity at selected positions in and around a target body (106) to provide improved heating of 

tude and phase weighting for the phased array transmit elements (104) of the hyperthermia applicator. A 915 MHz monopole 
^n/JZ^^ H app,lcato r for heati "8 brai » tuniors has an enclosed vessel including a plurality of monopole transmit 

tc TJ^^ZTJ, ™, 6 i 3S 3 C ' r ^ ar 3rC a 7 8y °" 3 g -T nd P,3ne Which has an aperture for Portioning the tumor Z proximity 
1? h! ™ P «° ! a 5! enna I el ^ ment ^- Adaptive focusing with non-invasive electric field probes is used to maximize the electric field 
at the tumor site. Parallel plate microwave waveguides are used to direct R.F. energy from the monopole phased array to the tu- 
mor site A microwave transmit and receive module generates amplitude and phase controlled transmit signals for exciting the 
monopole antenna elements, and receives passive microwave signals from the monopole antenna elements for taking " 
sive radiometry temperature measurements of the tumor site. B 
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SING and bmUJMfi HY°FR1 

MQNQPQI F PHASFP ARRAY AP P | frA-r^p 
COPYRIGHT 

Appendices A-D of the disclosure of this patent document contain 
5 materia, which Is subject to copyright protection. The copyright owner 
has no objection to the facsimile reproduction by anyone of the patent 
document or the patent disclosure, as it appears in the Patent and 
Trademark Office patent file or records, but otherwise reserves all 
copyright rights whatsoever. 
10 BACKGROUN D QFTHf [flyFJsjTjQN, 

The successful treatment of deep-seated malignant tumors within a 
patient is often a difficult task. The objective of the treatment is to reduce 
m size or compietely remove the tumor mass by one or more modalities 
available at the treatment facility. Common treatment modalities are 
1 5 surgery, chemotherapy, and x-ray therapy. One treatment modality used 
alone or in conjunction with one of the above modalities is "tissue 
heating", or hyperthermia. Hyperthermia can be considered as a form of 
h.gh fever localized within the body. A controlled thermal dose 
d.stribution is required for hyperthermia to have a therapeutic value 
20 Typical localized-hyperthermia temperatures required for therapeutic 
treatment of cancer are in the 43-45'C range. Norma, tissue should be 
kept at temperatures be.ow 43-C during the treatment. Typically 
hyperthermia is induced in the body by radio-frequency (RF) waves 
acoustic (u.trasound, waves, or a combination of both. One of the most 
25 difficult aspects of implementing hyperthermia, with either RF or 
ultrasound waves, is producing sufficient heating at depth 
Mu.tiple-applicator RF hyperthermia arrays are commonly used to provide a 
focused near-field main beam at the tumor position, Ideally, a focal region 
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should be concentrated at the tumor site with minimal energy delivered to 
surrounding normal tissue. 

In RF hyperthermia systems, the hyperthermia antenna beamwidth 
is proportional to the RF wavelength in the body. A small focal region 
5 suggests that the RF wavelength be as small as possible. However, due 
to propagation losses in tissue, the RF depth of penetration decreases with 
increasing transmit frequency. One of the major sjde-effects in heating a 
deep-seated tumor with a hyperthermia antenna is the formation of 
undesired "hot spots" in surrounding tissue. This additional undesired 
10 heating often produces pain, burns, and blistering in the patient, which 
requires terminating the treatment immediately. The patient does not 
receive anesthetics during the hyperthermia treatment in order to provide 
direct verbal feedback of any pain. Thus, techniques for reducing hot 
spots while maximizing energy delivered to the tumor site are desired in 
15 hyperthermia treatment. 

RF hyperthermia systems with electric field transmitting arrays, i.e., 
antenna arrays, in the frequency band of 60-2000 Mhz have been used to 
localize heating of malignant tumors within a target body. Phase control 
alone of the transmitting antennas of such an array has been used to 
20 synthesize therapeutic RF radiation patterns within a target body. 

Theoretical studies of adaptive control of individual antenna phase and 
power (transmit weights) has been used to maximize the tumor 
temperature (or RF power delivered to the tumor) while minimizing the 
surrounding tissue temperature (or RF power delivered to the surrounding 
25 tissue), invasive temperature measuring techniques have been used to 
optimize the radiation pattern within a target body. 

One commercially available hyperthermia annular phased-array 
antenna system is the Model BSD-2000, SIGMA-60 applicator, available 
from BSD Medical Corporation, Salt Lake City, Utah. This phased-array 
30 system fully surrounds the patient, placing the patient at the center of an 
annular array of dipole transmit antennas. By fully surrounding the patient 
with an annular phased-array, it is possible to obtain constructive 
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interference (o7 signal enhancement) deep within the target volume. This 
hyperthermia system uses a 60 cm array diameter with eight uniformly 
spaced dipole elements operating over the frequency band 60-120 MHz. 
The eight dipoles are fed as four active pairs of elements. There are four 
5 high-power amplifiers which drive the dipole pairs with up to 500 W 
average power per channel. Each of the four active channels has an 
electronically controlled variable-phase shifter for focusing the array. 
Temperature and electric-field probe sensors (both invasive and non- 
invasive) are used to monitor the treatment. A cool-water (5-40«»C) bolus 
10 between the patient and the phased-array is used to prevent excess 
heating of the skin surface. The water bolus is filled with circulating 
distilled water, which has a very low propagation loss. 

SUMMARY op tmp .fuv/pM-r^ N 
In accordance with the invention, adaptive nulling and/or focusing 
with non-invasive auxiliary probes is used to reduce or enhance the field 
Intensity at selected positions in and around the target body while 
maintaining a desired focus at a tumor thereby avoiding or reducing the 
occurrences of "hot spots" while enhancing heating of the tumor during 
ultrasonic or R.F. hyperthermia treatment. 
20 m general, in one aspect, the invention features a hyperthermia 

applicator having electric field radiators each coupled to a source of 
electric radiation through a controllable transmit weighting network to 
control the phase and amplitude of the electric field radiation transmitted 
by each radiator. The transmit weighting networks respond to feedback 
signals from a controller coupled to electric field probes which receive the 
electric field radiation from the radiators. The controller adjusts the 
feedback signals in response to the received electric field radiation so that 
the electric field radiation is minimized at the electric field probes. 

Preferred embodiments include a phased array of electric field 
radiators, and an annular array of electric field radiators for surrounding 
the target. 



25 



30 



SUBSTITUTE SHEET 



WO 93/00132 




PCT/US92/05464 




The electric field probes include probes placed non-invasively around the 
perimeter of the target where the electric field energy is to be minimized. 
In one embodiment, the target is modeled as an ellipse and the electric 
field probes are placed at the front, back, and on both sides of the ellipse. 



field probe, and the controller adjusts the feedback signals in response to 
the electric field radiation received by the secondary electric field probe so 
that the electric field radiation is maximized at the secondary probe. 
Embodiments include placing the secondary probe at the desired focus of 

10 the electric field radiation. 

In yet another aspect, the invention features the controller 
performing either a matrix inversion algorithm or a gradient search 
algorithm to adjust the feedback signals controlling the transmit weighting 
networks in response to the electric field energy received by the electric 

15 field probes. 

In general, in another aspect, the invention features a hyperthermia 
applicator for heating a target inside a body, having electric field radiators 
each coupled to a source of electric radiation through a controllable 
transmit weighting network to control the phase and amplitude of the 

20 electric field radiation transmitted by each radiator. The transmit 

weighting networks respond to feedback signals from a controller coupled 
to electric field probes placed outside the body which receive the electric 
field radiation from the radiators. The controller adjusts the feedback 
signals in response to the electric field radiation received outside the body 

25 so that the electric field radiation is controlled at the target inside the 
body. 

Preferred embodiments include a phased-array of electric field 
radiators, an annular array of electric field radiators for surrounding the 
target, and an array of monopole antenna elements for positioning nearby 
30 the target. 

In preferred embodiments of the monopole array, the monopole 
antenna elements are perpendicularly mounted to one side of an RF 



5 



In another aspect, the invention also features a secondary electric 
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reflecting groundplane. An RF reflecting screen is mounted perpendicular 
to the groundplane surface behind the monopole antenna elements to 
reflect RF energy from the monopole antenna elements toward the target. 
The ground plane includes an aperture for positioning the target on the 
5 same side of the ground plane as the monopole antenna elements. An 
enclosure surrounds the monopole antenna elements and provides a vessel 
for enclosing a bolus of fluid, such as deionized water, between the 
monopole antenna elements and the body. In other preferred 
embodiments another ground plane is provided above the monopole 
10 antenna elements to form a waveguide between the antenna elements and 

the target body. Further, multiple waveguides and monopole antenna 

arrays can be stacked. 

In other preferred embodiments of the monopole array, the 
monopole array antenna elements resonate between 800 and 1000 MHz 
1 5 and are arranged in a circular arc having a radius of between 5 and 20 cm. 
The body is a cranium and the target is a brain tumor. The radius of the 
monopole array circular arc is either the distance from the monopole 
antenna array to the center of the cranium, or to the target tumor. 

In still other preferred embodiments, the electric field probes are 
20 non-invasively placed along the perimeter of the body between the 

elements of the phased-array and the target. The controller adjusts the 
feedback signal, with a gradient search or matrix inversion algorithm, to 
minimize the difference In the electric field detected by adjacent electric 
field probes and thereby provide uniform electric field radiation into the 
25 body. 

In yet other preferred embodiments, the electric field probes are 
formed into an array non-invasively placed between the phased array and 
the target. The controller adjusts the feedback signal, with a gradient 
search or matrix inversion algorithm, to provide a particular electric field 
30 pattern across the electric field probe array and thereby focus radiation 
into the target. The electric field probe array elements are placed 
symmetrically with respect to a bisector line which runs from the target to 
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the phased array to bisect the phased array. The controller adjusts the 
feedback signal to balance the electric field pattern with respect to the 
bisector line, and to minimize the difference in the electric field detected 
along the bisector line. 
5 Thus, the present invention offers the advantages of allowing 

effective hyperthermia treatment to be applied to deep-seated tumors 
within the body while reducing or eliminating hot-spot formation on the 
surface of the body which interferes with the treatment. Another 
advantage is that hot spots are eliminated quickly by sensing and adjusting 
10 the E-field radiation in the vicinity of the expected hot spot rather than by 
measuring the temperature rise of the tissue after heating has already 
occurred. Still another advantage is that the E-field sensing probes may be 
located on the surface of the target rather than having to be invasively 
placed within the target body. A further advantage is that the E-field 
15 radiation can be focused on a target inside a body using E-field sensing 
probes non-invasively placed outside the body to maximize heating of the 
target tissue. 

BRIEF DESCRI PTION OF THE DRAWING 
The foregoing and other objects, features and advantages of the 
20 invention will be apparent from the following more particular description of 
preferred embodiments of the invention, as illustrated in the accompanying 
drawings in which like reference characters refer to the same parts 
throughout different views. The drawings are not necessarily to scale, 
emphasis instead being placed upon illustrating the principles of the 
25 invention. 

FIG. 1 is a perspective view of an RF annular array hyperthermia 
system featuring the adaptive nulling of this invention. 

FIG. 2 is a cross-sectional representation of the annular array of FIG. 

1. 

30 FIG. 3 is an analytical model of the cross-sectional representation of 

FIG. 2. 
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FIG. 4 is 8 simulated thermal profile of the analytical model of FIG. 3 
without the adaptive nulling of this invention. 

FIG. 5 is a simulated thermal profile of the analytical model of FIG. 3 
with the adaptive nulling of this invention. 
5 FIG. 6 is a schematic diagram of the adaptive hyperthermia array 

and array controller of FIG. 1 . 

FIG. 7 is a schematic diagram of an analytical model of an adaptive 
array for simulating the hyperthermia array of FIG. 1 . 

FIG. 8 is a block diagram detailing the sample matrix inversion 
10 algorithm derived from the adaptive hyperthermia array model of FIG. 7. 

FIG. 9 is a block diagram of the sample matrix inversion algorithm 
performed by the hyperthermia array controller of FIG. 6. 

FIG. 10 is a scatter diagram of transmit weights used in deriving the 
gradient search adaptive hyperthermia algorithm. 
1 5 FIG. 1 1 Is a diagram showing the derivation of the gradient search 

directions. 

FIG. 12 is a block diagram of the gradient search performed by the 
hyperthermia array controller of FIG. 6. 

FIG. 13 is a schematic diagram of the analytical model of FIG. 7 
20 redrawn to simplify derivation of method of moments analysis. 

FIG. 14 is a schematic diagram of an equivalent circuit model for 
simulating an auxiliary probe. 

FIG. 15 is a schematic diagram of a thermal conductivity model for 
simulating hyperthermia heating within a target. 
25 FIG. 1 6(a) is a block diagram detailing the simulation model of the 

hyperthermia array of FIG. 1. 

FIG. 16(b) is a table of values used in the simulation model of FIG 

16(a). 

FIG. 17 is a schematic diagram of the transmit antenna array and 
30 auxiliary probe array geometries for the simulation model of FIG. 1 6(a). 
FIG. 1 8 is a diagram of the simulated E-field for the simulation 
model of FIG. 16(a) prior to adaptive nulling. 
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FIG. 19 is a diagram of the simulated profile in 1 dB steps for the E- 
fieldof FIG. 18. 

FIGS. 20 and 21 are diagrams of the simulated E-field profile of FIG. 
19 taken along the x- and z- axes, respectively. 
5 FIG. 22 is a diagram of the simulated E-field for the simulation 

model of FIG. 1 6(a) after adaptive nulling. 

FIGS. 23 and 24 are diagrams of the simulated E-field of FIGS. 18 
and 22 taken along the x- and z- axes, respectively. 

FIG. 25(a) is a diagram of the simulated E-field before and after 
1 0 adaptive nulling taken longitudinally in the y direction along the x = 1 5, 
z=0 cm line of the geometry shown in FIG. 25(b). 

FIG. 26(a) is a diagram of the simulated E-field before and after 
adaptive nulling taken longitudinally along the y axis (x=0, z=0 cm) of 
the geometry shown in FIG. 26(b). 
15 FIGS. 27(a) and 27(b) are graphs showing the transmit weight 

amplitude and phase, respectively, before and after adaptive nulling. 

FIG. 28 is a graph showing the channel correlation matrix 
eigenvalues. 

FIG. 29 is a diagram of the simulated target temperature profile for 
20 the E-field of FIG. 18 prior to adaptive nulling. 

FIGS. 30 and 31 are diagrams of the temperature profile of FIG. 29 
taken along the x- and z- axes, respectively. 

FIG. 32 is a diagram of the simulated target temperature profile for 
the E-field of FIG. 22 after adaptive nulling. 
25 FIGS. 33 and 34 are diagrams of the temperature profile of FIGS. 

29 and 32 taken along the x- and z- axes, respectively. 

FIG. 35 is a diagram of the simulated target temperature profile for 
the E-field of FIG. 18 prior to adaptive nulling. 

FIG. 36 is a diagram of the simulated target temperature profile for 
30 the E-field of FIG. 22 after adaptive nulling. 

FIGS. 37 and 38 are diagrams of the temperature profile of FIG. 36 
taken along the x- and z- axis, respectively. 
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FIG. 39 is a cross-sectional view of an annular phased array 
hyperthermia system and saline phantom used for gathering experimental 
adaptive nulling and focusing data. 

FIG. 40 is a diagram of the E-field amplitude measured at a single 
5 null site, the tumor site, and a reference site, versus gradient search 
iteration for the experimental hyperthermia system and phantom of FIG 
39. 

FIG. 41 is a diagram of the simulated E-field power versus gradient 
search iteration at the single null site of FIG. 40. 
10 FIG. 42 is a cross-sectional view of a beef phantom used in place of 

the saline phantom in the experimental system of FIG. 39. 

FIG. 43 is a diagram of the E-field amplitude measured at a single 
null site, the tumor site, and a reference site, versus gradient search 
iteration for the experimental hyperthermia system of FIG. 39 using the 
1 5 beef phantom of FIG. 42. 

FIG. 44 is a diagram of the measured temperature versus time at 
the beef phantom tumor site and the single null site of FIG. 42. 

FIG. 45 is a diagram of the E-field amplitude measured at two null 
sites and the tumor site, versus gradient search iteration for the 
20 experimental hyperthermia system and phantom of FIG. 39. 

FIG. 46 is a diagram of the E-field amplitude measured at a focus on 
the surface of the saline phantom versus gradient search iteration for the 
experimental hyperthermia system and phantom of FIG. 39. 

FIG. 47 is a perspective view of an RF monopole array hyperthermia 
25 system for treating brain tumors, featuring the adaptive focusing of this 
invention. 

FIG. 48 is a top-view of the monopole phased array hyperthermia 
system of FIG. 47 

FIG. 49 is a cross-sectional side view of the RF monopole phased 
30 array hyperthermia system of FIG. 47. 

FIG. 50 is a schematic diagram of the monopole phased array 
hyperthermia system of FIG. 47 and an array controller. 
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FIG. 51 is a schematic diagram of the monopole phased array 
geometries for the hyperthermia system of FIG. 47. 

FIG. 52 is a schematic diagram of the monopole phased array and 
electric field probe array geometries for the hyperthermia system of FIG. 



47. 



FIG. 53 is a diagram of a simulated E-field pattern for the monopole 
phased array geometry of FIG. 52. 

FIG. 54 is a diagram of the simulated E-field pattern of FIG. 53 
taken parallel to the x- axis along the line 2 = 5.08 cm. 
10 FIG. 55 is a diagram of the simulated temperature profile for the 

simulated E-field pattern of FIG. 53. 

FIG. 56 is a schematic diagram of the monopole phased array " 
hyperthermia system of FIG. 47 with an electric field probe array for 
generating uniform RF illumination over a large area of a target. 
1 5 FIG. 57 show the monopole phased array system of FIG 48 

including the addition of a top ground plane surface forming a parallel plate 
waveguide. 

FIGS. 58A shows a cross-sectional view of the wave guide of FIG. 



57. 



20 



25 



FIG. 58B shows a cross-sectional view of an alternative 
embodiment of the waveguide of FIG. 58A having diverging surfaces. 

FIG. 58C shows a cross-sectional view of another alternative 
embodiment of the waveguide of FIG. 58A having flared surfaces forming 
a horn. H 

FIG. 59 shows a cross-sectional view of another alternative 
embodiment of the waveguide of FIG. 58A having multiple stacked 
monopole phased array antennas and associated stacked waveguides. 

FIG. 60 is a block diagram of a microwave transmit and receive 
module for use with the monopole phased array system of FIG. 47. 
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PETAIIEP PFftCRIPTIQN OF THF PREFERRED mSQQIMEKL 
Apparatus 

Referring to FIG. 1. there is shown a hyperthermia annular phased- 
array system 100 having improved "hot spot" characteristics achieved by 
5 utilizing the focused near-field adaptive nulling apparatus of this Invention. 
An annular hyperthermia phased-array applicator 102, energized by a 
hyperthermia array controller 101. has a plurality of dipole transmit 
antenna elements 104 placed around a patient to be treated, or target 
body 106. The dipole antenna elements are uniformly disposed around the 
10 patient. Each dipole antenna element is oriented parallel to the other 

dipole antennas and parallel to a longitudinal axis A-A passing through the 
center of a cylinder defined by applicator 102. The patient is positioned 
within the hyperthermia phased array applicator 102 such that the deep- 
seated tumor to be treated 107 is at the approximate center, or focus, of 
15 the phased array applicator. A water-bolus 105 is provided between the 
patient and the phased array applicator to control the temperature of the 
patient's skin. Phased-array applicator 102 therapeutically illuminates the 
target body 106 with electric field (E-field) energy radiated by dipole 
antenna elements 104 focused on tumor 107 deep within the body. 
20 An example of a deep-seated tumor is cancer of the prostate. The 

tumor volume often has a decreased blood flow which aids in heating the 
tumor, compared to normal tissue for which heat is carried away by 
normal blood flow. In practice, undesired high-temperature regions away 
from the focus can also occur on the skin and inside the volume of the 
25 target body. For example, scar tissue, which has a decreased blood flow 
rate, will tend to heat up more rapidly than normal tissue having normal 
blood flow. 

In the adaptive hyperthermia array of this invention, electric-field 
nulls are used to reduce the power delivered to potential hot spots. 
30 Computer simulations, described herein, establish that non-invasive field 
probes, or sensors, 112 placed on the surface of the target can be used to 
eliminate hot spots interior to the target tissue. With the adaptive 
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hyperthermia phased-array described herein, RF energy nulls are adaptively 
formed to reduce the electric field energy delivered to these potential hot 
spots. As will be shown, the energy nulls achieved by the adaptive nulling 
apparatus of this invention are both invasive to the target, i.e., extend into 
5 the target body, and non-invasive to the target, i.e., on the surface of the 
target. 

Referring to FIG. 2, there is shown a schematic cross-sectional 
representation of an embodiment of an eight-element hyperthermia 
phased-array applicator 102 of FIG. 1. Phased-array applicator 102 has 
10 transmit antennas 104, through 104 8 , arranged symmetrically surrounding 
a human body target 106' at the prostate level. 

An analytical model of the embodiment of FIG. 2 is shown in FIG. 3. 
Here, an elliptical phantom target 106 is used to model the prostate-level 
cross section of the human body 106'. The center 107 of the elliptical 
15 phantom models the location of the prostate tumor to receive 

hyperthermia treatment, i.e., the focus of RF energy for the phased array 
applicator 102. Water bolus 105 is assumed to surround the target body 
106, and is treated as a homogeneous medium for analysis purposes. 

Four auxiliary RF E-field probes, or sensors, 112, through 112 4 , i.e., 
20 receiving antennas, are placed around the perimeter of the target to model 
non-invasive probes placed on the skin of the human body target. Each 
auxiliary probe 112, through 1 12 4 has a corresponding null zone 120, 
through 1204, respectively, centered at each auxiliary probe and extending 
into the elliptical target region 106. Each null zone indicates an area in 
25 which undesired "hot spots" are reduced or eliminated. The width of each 
null zone is directly related to the strength of each null. The strength of 
each null (sometimes referred to as the amount of cancellation) is directly 
related to the signal-to-noise ratio at the probe position (SNR P ). A low 
SNR P indicates a large amount of nulling (strong null), and a high SNR p 
30 indicates a small amount of nulling (weak null). The resolution, or 

minimum spacing, between the focus 107 and any null position is normally 
equal to the half-power beamwidth of the transmit antenna. Resolution 
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may be enhanced somewhat by using weak nulls whenever the separation 
between the null and focus is closer than the half-power beamwldth. 

Referring to FIG. 4, there is shown the results of a simulation of 
the thermal distribution inside the target body 106 for the hyperthermia 
5 ring array applicator 102 of the analytical model of FIG. 3, without 
adaptive nulling, transmitting into the target body. For simulation 
purposes, target body 106 is assumed to be a homogeneous elliptical 
region, and the RF energy from the array is focused at the center of the 
ellipse 107, simulating the tumor site. No adaptive nulling is used. The 
10 contour lines of the thermal distribution represent isotherms having the 
indicated temperature in degrees Celsius (°C), and are spaced at 2°C 
intervals. The simulation shows that the focus is expected to have a 
temperature of approximately 46°C, while two undesired "hot spots" 122 
and 1 24 to the left and right of the focus, respectively, are expected to 
15 have temperatures of approximately 42 °C. 

FIG. 5 shows a simulated thermal distribution for the model of FIG. 
3 where the adaptive nulling methods of this invention are applied. 
Comparison of FIG. 5 with FIG. 4 show that the "hot spots" 122 and 124 
are essentially eliminated, no new "hot spots" have been produced within 
20 the target body, and the peak temperature induced at the focus is still 
approximately 46°C. 

Referring to FIG. 6, a generalized schematic of the non-invasive 
adaptive-nulling hyperthermia system of FIG. 3 includes hyperthermia 
transmitting phased array applicator 102 having a plurality of transmitting 

25 antenna elements 1 04 n , where n « 1 N, surrounding target body 1 06 for 

focusing RF energy at focus 107 within the target body. Phased array 
applicator 1 02 is energized by an RF energy source 1 08 which is 
distributed to and drives each transmit antenna element 104 n through a 
corresponding transmit weighting function 1 10 n , each having a 

30 corresponding weight w n . Each weighting function w n may affect the gain 
and phase of the RF energy fed to its corresponding antenna 104 n in the 
array, i.e., w n represents a complex weighting function. Each weighting 
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function 1 10 n may be implemented by a voitage controlled RF amplifier 
and a voltage controlled RF phase shifter. An amplitude control voltage 
representing the amplitude component of transmit weight w n is fed to the 
voltage controlled amplifier, and a phase control voltage representing the 
Phase of transmit weight w n is fed to the voltage controlled phase shifter. 
Target body 106 has a plurality of E-field auxiliary probes 1 12 m , 

where m = 1 N aux , i.e., receiving antennas, positioned at various ^ 

locations on the surface of the body for sampling the E-field at each 
particular location. Another receiving probe 1 15 may be placed at the 
1 0 desired focus 1 07 of the array. 

Receiving probes 1 12 m and 1 15 each drive an input to an RF 
receiver 114. The transmit amplitude and phase weights of each 
weighting function w n are fed to the receiver 1 14 through lines 103 n and 
are used to find the transmit level of each transmit element 104 The 
15 outputs of receiver 1 14 represent the auxiliary probe-received complex 
voltages v u v 2r . . the focus pro5e . rece|ved comp|ex 

and the transmit level of the phased array. The receiver outputs drive the 
•nputs of a signal processor 116, which applies a nulling algorithm to 
adjust the weighting functions w n and thereby null, or minimize, the RF 
s,gna, received by each receiving probe 1 12 m , i.e., minimize the SNR P at 
each probe. p 

To generate the desired field distribution in a clinical adaptive 
hyperthermia system, the receiving probes are positioned as close as 
possible to the focus (tumor site, and to where high temperatures are to 

25 be avoided (such as near the spina, cord and scar tissue,. For an annular 
array configuration the receiving probes can be located non-invasively on 
the surface (skin, of the target, .nitial.y, the hyperthermia array is focused 
to produce the required field intensity at the tumor. An invasive probe 
may be used to achieve the optimum focus at depth. To avoid undesired 

30 hot spots, It is necessary to minimize the power received at the desired 
null positions and to constrain the array transmit weights w n to deliver a 
requ.red amount of transmitted or focal region power. 
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Signal processor 116 performs either a sample matrix inversion 
(SMI) algorithm or a gradient search algorithm on the signals output from 
receiver 114 and updates the adaptive array weights w n (with gain g and 
phase 0} to rapidly (within seconds) form the nulls at the auxiliary probes 
5 before a significant amount of target heating takes place. With this 
adaptive system, it is possible to avoid unintentional hot spots in the 
proximity of the auxiliary probes and maintain a therapeutic thermal dose 
distribution at the focus (tumor). 

Signal processor 1 16 may also perform a maximizing algorithm to 
10 maximize energy at the focus 1 07. The focus probe 1 1 5 is invasively 
Placed at the desired focus 107, and used to generate a maximum signal, 
or signal-to-noise ratio (SNR F ), at the tumor site. RF receiver 1 14 makes 
an amplitude and phase measurement on the output signal from invasive 
probe 1 15 for each transmit antenna element 104 n radiating one at a time. 
15 Signal processor 1 1 6 processes these measurements and feeds back 
weight command signals to the transmit weighting functions 1 10 n to 
calibrate or phase align the transmit channels to thereby maximize" the 
SNR F . or RF power, at the invasive focal point probe. If receiver 1 14 
makes amplitude-only measurements from invasive focus probe 115, then 
20 a gradient search technique may be applied by the signal processor with all 
elements transmitting simultaneously to maximize the SNR F at the invasive 
focal point probe. 

Theoretical Formulatinn ftf Nullinn Alpnrft^ 
FIG. 7 shows an analytical model of a hyperthermia phased-array 
25 antenna system 200, paralleling the generalized hyperthermia phased-array 
antenna system 100 of FIG. 6, Illustrating the principles of the near-field 
adaptive nulling technique of this invention. The phased-array antenna 
system 200 includes a hyperthermia transmitting antenna array 202 

having a plurality of transmitting antennas 204 n> where n = 1 N for 

30 focusing RF energy at a desired focus 207 in the near field of the antenna. 
Antenna array 202 is energized by an RF energy source 208 which drives 
a power divider 209. Power divider 209 has one output for driving each 
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antenna 204 n a corresponding transmit weighting function 210 n , each 
having a corresponding transmit weight w n . It is assumed here that each 
weighting function w n may affect the phase of the RF energy fed to its 
corresponding antenna 204 n in the array. A calibration E-field probe 212, 
5 or focus probe antenna, is positioned at focus 207 for sampling the E-field 
at that location. 

It is assumed that the hyperthermia phased-array antenna 200 is 
focused (as it normally is) in the near field and that a main beam 220 and 
possibly sidelobes 222 are formed in the target. In general, phase and 
10 amplitude focusing is possible. It is assumed that phase focusing alone is 
used to produce the desired quiescent main beam, i.e., weighting 
functions w n affects only the phase of the RF signal driving each antenna. 
The signal received by the calibration probe can be maximized by adjusting 
the phase weighting functions w n so that the observed transmit antenna 
1 5 element-to-element phase variation is removed, i.e., all transmit antennas 
appear to be in-phase when observed from the focus. 

One way to achieve phase coherence at the focus in a numerical 
simulation is to choose a reference path length as the distance from the 
focus to the phase center 224 of the array. This distance is denoted r F 
20 and the distance from the focus to the nth array transmit antenna element 
is denoted r/. The voltage received at the calibration probe 212 (located 
at focus 207) due to the mh array element may be computed using the 
-method of moments", as described below. To maximize the received 
voltage at the calibration probe output, it is necessary to apply the phase 
25 conjugate of the signal observed at the calibration probe, due to each 
array transmit antenna element, to the corresponding element at the 
transmit array. The resulting near-field radiation pattern will have a main 
beam and sidelobes. The main beam will be pointed at the array focal 
point, and sidelobes will exist at angles away from the main beam. 
30 Auxiliary probes can then be placed at the desired null positions in the 
quiescent sidelobe region. These sidelobes occur where tissue hot spots 
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are Hkely to occur, and they are nulled by one of the adaptive nulling 
algorithms described below. 

Adaptive Transmit Arra y Formulation 

Considering again the hyperthermia array and probe geometry 
5 shown in FIG. 7, the hyperthermia transmit antenna array 202 typically 
contains N identical transmit antenna elements 204. The number of 
adaptive channels is denoted M, and for a fully adaptive array M =/V. The 
ideal transmit weights w n (a complex voltage gain vector) are assumed in 
the computer simulation, with w = (w,,^...,^ denoting the adaptive 
10 channel weight vector as shown in FIG. 6. (Superscript T means 

transpose). To generate adaptive nulls, the transmit weights (phase and 
gain) are controlled by either the Sample Matrix Inversion (SMI) algorithm 
or a gradient search algorithm. The SMI algorithm has the flexibility to 
operate in either open- or closed-loop feedback modes; the gradient search 
15 algorithm operates only in a feedback mode. 

Sample Matrix Inversion (fl M M Alnorlth™ 

For the SMI algorithm, the fundamental quantities required to fully 
characterize the incident field for adaptive nulling purposes are the 
adaptive channel cross correlations. To implement this algorithm it is 
20 necessary to know the complex received voltage at each of the auxiliary 
probes. For example, the moment-method formulation (described below) 
allows computation of complex-received voltage at each of the auxiliary 
probes. 

FIG. 8 is a block diagram showing the SMI algorithm applied to the 
25 adaptive hyperthermia phased-array of FIG. 7. and the derivation of 
performance measures to quantify computer simulation results. Four 
performance measures are used to quantify the computer simulations: 
electric-field distribution E(x,y. Z ), channel correlation matrix eigenvalues 
A * k = 1 N ' 8da P tiv e transmit weights w., and interference cancellation 
30 C. The calculation of these performance measures is described in detail 
below. 
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Assuming a spherical wavefront is incident at an Ah probe antenna 
226 due to each of the N array transmit antenna elements 204 n (radiating 
one at a time with a unity-amplitude reference signal), the result is a set of 
probe-received complex voltages denoted vf,v£ • • .,vj after a gain 
adjustment 250. The cross correlation of the received voltages due to 
the /mh and nth transmit antenna (adaptive transmit channel) at the Ah 
probe is given by 



10 



RL=my m vd . (i) 

where * means complex conjugate and E(.) means mathematical 
expectation. (Note: for convenience, in Equation (1) the superscript / in v m 
and in v n has been omitted.) Because v m and v„ represent voltages of the" 
same waveform but at different times, R m ' n is also referred to as an 
autocorrelation function. 

In the frequency domain, assuming the transmit waveform has a 
band-limited white noise power spectral density (as commonly assumed in 
1 5 radar system analysis), Equation (1 ) can be expressed as the frequency 
average 



*L=yfaf)v;(f)df , (2 , 

where B = f 2 -f f is the nulling bandwidth, or bandwidth of frequencies 
applied by the hyperthermia treatment, and f is the transmit frequency of 
the hyperthermia array. It should be noted that vjf) takes into account 
20 the transmit wavefront shape, which is spherical for the hyperthermia 
application. For the special case of a continuous wave (CW) transmit 
waveform, as normally used in hyperthermia, the cross correlation reduces 
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where f Q is the transmit frequency of the hyperthermia array. 

Next, the channel correlation matrix, or interference covariance 
matrix, denoted R is determined 252. (Note: in hyperthermia, interference 
is used to refer to the signals received at the auxiliary probes. The 
undesired "hot spots" can be thought of as interfering with the therapy.) 
If there are N gux independent desired null positions or auxiliary probes, the 
/Wprobe channel correlation matrix is the sum of the channel correlation 
matrices observed at the individual probes. That is, 



10 



15 



where F, is the sample channel correlation matrix observed at the /th probe 
and / is the identity matrix used to represent the thermal noise level of the 
receiver for simulation purposes. 

Prior to generating an adaptive null, the adaptive channel weight 
vector, w, is chosen to synthesize a desired quiescent radiation pattern 
When nulling is desired, the optimum set of transmit weights to form an 
adaptive null (or nulls), denoted w m , is computed 254 by 



(5) 



where means inverse and w q is the quiescent weight vector. During 
array calibration, the normalized quiescent transmit weight vector, with 
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transmit element 204, radiating, is chosen to be w q =n.O,0,-- .,0) T , i.e., 
the transmit channel weight of element 204, is unity and the remaining " 
transmit channel weights are zero. Similar weight settings are used to 
calibrate the remaining transmit elements. For a fully adaptive annular 
5 array focused at the origin in homogeneous tissue, the normalized 
quiescent weight vector is simply w q =f1. /,/,•• fjT Commonly, the 
weight vector is constrained to deliver a required amount of power to the 
hyperthermia array or to the tumor. For simplicity in the computer 
simulation used to analyze the hyperthermia array, the weights are 
tO constrained such that 



£ KM . (6) 

where w„ is the transmit weight for the mh element. It should be noted 

that in the computer simulations, the electric field due to the normalized 

weight vector is scaled appropriately to deliver the required amount of 

power to the tissue so that a desired focal-region temperature level is 

15 achieved after t minutes. The summation of power received at the probes 
is given by 



P=**Rw , (7) 



where * means complex conjugate transpose. The 
signal-plus-noise-to-noise ratio for the auxiliary probe array, denoted SNR^ 
20 « computed as the ratio of the auxiliary probe array output power (defined 
.n Equation (7)) with the transmit signal present, to the probe array output 
power with only receiver noise present, that is. 
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Next, the adaptive array cancellation ratio indicative of the null strength, 
denoted C, is determined 255. C is defined here as the ratio of the 
summation of probe-received power after adaptation to the summation of 
probe-received power before adaptation (quiescent); that is, 

C= — • (9) 



A large amount of cancellation indicated by a large value for C indicates a 
strong null, while a small amount of cancellation indicated by a small value 
for C indicates a weak null. Substituting Equation (7) into Equation (9) 
yields 



^ ' (1 °> 



10 Next, the channel correlation matrix defined by the elements in 

Equations (2) or (3) is Hermitian (that is, R=r\ which, by the spectral 
theorem, can be decomposed 256 in eigenspace as 



u 

(id 



where A k , * = 1,2, • • >,M are the eigenvalues of R, and k= 1,2, • . 
;M are the associated eigenvectors of R. The channel correlation matrix 
15 eigenvalues M,,^ • • :A M ) are a convenient quantitative measure of the 
use of the adaptive array degrees of freedom. The amplitude spread 
between the largest and smallest eigenvalues is a quantitative measure of 
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the dynamic range of the interference (hot spot) signals. FIG. 9 is a block 
diagram of the sample matrix inversion algorithm implemented by the 
signal processor 1 1 6 of FIG. 6. Receiver 1 14 generates probe-received 
complex voltage vector vj, v£ -vi^ for the hh auxiliary probe. The 
5 signal processor generates 280 the transmit channel correlations R^ 
defined by equation (3), and sums 282 them to form the channel 
correlation matrix R defined by equation (4). Next, the signal processor 
multiplies 284 the inverse of the channel correlation matrix fT 1 by the 
quiescent transmit weight vector w g to form the new adapted transmit 
10 weight vector w s containing the adapted transmit weights fed back to the 
transmit weight networks 1 10 n of FIG. 6. 

Gradient Search Alflflrifhrn 

Under conditions where only the probe received voltage amplitude is 
measured, it is appropriate to consider a gradient search algorithm to 
15 minimize the interference power at selected positions. The gradient search 
is used to control the transmit weights w„ iteratively such that the RF 
signal received by the probe array is minimized. The transmit array 
weights (gain and phase) are adaptively changed in small increments and 
the probe array output power is monitored to determine weight settings 
20 that reduce the output power most rapidly to a null. 

Consider J sets of N transmit weights that are applied to adaptive 
hyperthermia phased array applicator 102 of FIG. 6. In terms of adaptive 
nulling, the optimum transmit weight settings (from the collection of J sets 
of N transmit weights) occur when the SNR P is minimized. Equivalently, 
25 the total interference power received by the auxiliary probe array, denoted 
p^ c , is to be minimized. For notational convenience let a figure of merit F 
denote either the SNR P orpP 0 and employ a gradient search to find the 
optimum transmit weights to minimize F, that is, 



^min(A5/*1,2,...,J . 



(12) 
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The transmit weight settings for which F opt occurs yields the closest 
approximation to the optimal transmit weights determined by using the 
sample matrix inverse approach described above. 

FIG. 10 shows an amplitude and phase scatter diagram for the N 
complex transmit weights w n at the/th configuration. I.e., the /th set of 
weights tried. The nth transmit weight in the/th configuration of transmit 
weights is denoted 



where A nJ is the transmit weight amplitude distributed over the ranoe A ■ 
to A max and <*>„/ is the transmit weight phase distributed over the range -n 
10 to n radians. 

Referring also to FIG. 11, it is desired to find the values of 
amplitude and phase for each of the N transmit weights such that the 
figure of merit F (SNR p or p™) is minimized. When the figure of merit is 
minimized, adaptive radiation pattern nulls will be formed at the auxiliary 
15 probe positions. 

Assuming an initial setting of the N transmit weights such as those 
selected to focus the radiation pattern on a tumor, the weights are 
adjusted by dithering them until the optimum figure of merit is achieved. 
It is desired to find the collective search directions for the N transmit 

20 weights such that /-decreases most rapidly. That is, weights are selected 
so that the directional derivative is minimized at iA,*,), where Aj and O y 
are vectors representing the transmit amplitude weights and transmit 
phase weights, respectively, for the /th configuration. 

The directional derivative of Fj is expressed in terms of the 

25 amplitude and phase changes of the transmit weights as 
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where 3 means partial derivative, and r An/ are the {A, <P) directions for 
which Fj is decreasing most rapidly. The directions r^ n] are 
constrained by 



r<cKM • (15) 



It is desired to minimize DfFjf subject to the above constraint equation. 
5 Using Lagrange multipliers it is possible to construct the Lagrangian 

function 



where G is a constant to be determined. The requirement that L. be an 
extremum implies 



a^ = a^- 2G/ ^=o (17) 

and 

BL dF. 

a^a*^ 2 ^ 0 (18) 

10 or that 
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and 



(19) 



(20) 



Squaring equations (19) and (20) and invoking equation (15) yields 



(21) 



thus. 



2 N 



n-1 



3V 



(22) 
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Substituting this expression for G in equations (19) and (20) gives 



v- 



dFj 



n-1 dAj BQj 



(23) 



and 




The minus sign was chosen corresponding to the direction of maximum 
function decrease. This choice of minus sign in equation (22) enforces 
5 nulls in the hyperthermia array radiation pattern. Alternatively, if the 
positive sign in equation (22) is selected, then the gradient directions can 
be used to maximize the figure of merit for the purposes of focusing at an 
invasive probe at the tumor site, i.e., maximize the SNR F . This may be 
used, for example, to determine the quiescent transmit weight vector w . 
1 0 Thus, two gradient searches may be performed to optimize the radiation" 
pattern of the hyperthermia array. The first to produce a peak or focused 
radiation pattern at the tumor, and the second to form the desired nulls at 
the auxiliary probes. Furthermore, these two gradient searches may be 
implemented as a single, combined gradient search constrained to 
1 5 maximize the radiation pattern at the focus and minimize the radiation 

pattern at the desired nulls. The combined gradient search is implemented 
by minimizing the figure of merit defined as the ratio of the power received 
at the auxiliary probes to the power received by the probe at the focus. 
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The partial derivatives 




dF. dF. 

~d*j "^"'N (25) 

represent the gradient directions for maximum function decrease. Since 
the figure of merit F cannot be expressed here in analytical form, the 
partial derivatives are numerically evaluated by using finite differences. 
5 Thus, we write 



dF, A^ 



and 



a*„ 2a*^ (27) 



are 



where as shown in FIG. 2 the figure of merit differences 

AV/V^V^/V^^ (28) 

and 

A/ W , '3fW A »^-'5CW A »^ (29) 

and LA nJ ar\6 Acp^are assumed to be small increments. 
10 We will assume that the increments LA nJ and A<p„. are independent of the 
configuration number and element number, that is, 
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m 



AA^AA 



(30) 



and 



(31) 



Substituting equations (26), (27), (30) and (31) in equations (23) and (24) 
gives the desired result for the search directions 



5 and 



LA 



n-1 



AA A* 



(32) 



AF. 



A4> 



N 



n-1 



AA A* 



(33) 



Equations (32) and (33) are used to compute the new amplitude and phase 
settings of the flf+ r;th transmit weight configuration according to 
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A^^A^AAr^ 



(34) 



and 



(35) 



25 



In practice, ft may be necessary to keep one of the transmit weights fixed 
(in amplitude and in phase) during the gradient search to guarantee 
convergence. 

5 FIG. 1 2 is a block diagram of the gradient search algorithm 

implemented by the signal processor 1 16 of FIG. 6. Each of the N 
transmit antennas 104 n of phased array applicator 102 (FIG. 6) is driven 
through its corresponding weighting network 110 n which applies complex 
transmit weights w nj at theyth configuration of the weights. The transmit 
10 antennas induce a voltage across the fth probe antenna 1 12, at the 

corresponding input to receiver 1 14 (FIG. 6). Receiver 114 amplifies the 
signa. received from the fth probe by gain a,, to produce voltage amplitude 
vector } vft, \ v£ u . . ., | Vf ].x at the receiver output 

The voltage amplitude vector is input to signal processor 116 which 
performs the gradient search. For any initial configuration (/= 1 , of the 
transmit weights w nj , the signal processor causes each weight to be 
dithered by a small amount in amplitude, LA n] , and phase, A<J>„.. Each 
transmit weight is dithered Independent of the other transmit weights 
which remain in their /th configuration state. Received voltage vectors 
Wi'Mvjl- • -,1^1. i.e., are stored and used to calculate the resulting 
figure of merit F nJ 300 for each dithered condition, the figure of merit 
being the power received by the auxiliary prove array. The figure of merit 
•s a rectangular matrix of dimension N x 4, where the dimensionality of 
four is due to the plus and minus dithering of both of the amplitude and 
phase. The figure of merit differences AF AnJ and AF 0nJ caused by 



15 



20 
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dithering the amplitude and phase, respectively, are calculated according 
to equations (28) and (29). The gradient search directions r AnJ and r 9nj , 
based upon minimizing the auxiliary probe array received power, are then 
determined 302 from the figure of merit differences according to equations 
(32) and (33), respectively. The resulting search directions are used to 
update 304 transmit weights w n] to the (j+ 1)xb configuration transmit 
weights w n ff+ v according to equations (34) and (35). The transmit 
weights w n ff+v are sent to update the transmit weighting networks 1 10 n , 
and the process is repeated. The final adaptive weight vector w a is 
achieved when the ff+ J)th transmit weight configuration has converged. 
Convergence is expected to occur within several hundred iterations 
depending on the dither step size AA and A(p. 

It is understood that other forms of gradient searches exist which 
can be used to update the transmit weights toward convergence. Another 
15 such gradient search approach, where the step sizes LA and A<t> are 

computed at each iteration, is described by D.J. Farina and R.P. Flam, "A 
Self-normalizing Gradient Search Adaptive Array Algorithm", IEEE 
Transactions on Aerospace and Electronic Systems, November 1991, Vol. 
27, No. 6, pp 901-905. 



20 



COMPUTER SIMUL ATION O F A D APTIVE Mill I IM G HYPFRTHFPMi a 
Moment-Method Fnrmi ilatinn 

Referring again to FIG. 8, a method of moments formulation 258 is 
used to compute the probe-received voltages in Equation (2) due to the 
transmitting hyperthermia phased-array antenna in an infinite 

25 homogeneous conducting medium. The medium is described by the three 
parameters », € , and a. which are discussed below. The formulation given 
here is analogous to that developed under array-receiving conditions for an 
adaptive radar. The software used to analyze a hyperthermia array is 
based on the receive-array analogy but the theory presented below is 

30 given in the context of a transmit array. 
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An antenna analysis code (WIRES) originally developed by J.H. 
Richmond is capable of analyzing antenna or radar cross section problems. 
See., J.H. Richmond, "Computer program for thin-wire structures In a 
homogeneous conducting medium-. Ohio State University, ElectroScience 
5 Laboratory, Technical Report 2902-12, August 1973; and, J.H. Richmond, 
-Radiation and scattering by thin-wire structures in a homogeneous 
conducting medium (computer program description)-, IEEE Trans. 
Antennas Propagation, Vol. AP-22, no. 2, p.365, March 1974. WIRES 
was modified to analyzing the near-field and far-field adaptive nulling 
10 Performance of thin-wire phased arrays in free space. A new version of 
the thin-wire code that can analyze adaptive hyperthermia arrays in an 
infinite homogeneous conducting medium was written to conduct the 
adaptive hyperthermia simulation discussed below. The new version of 
the thin-wire code is attached as Appendix A. 
1 5 WIRES is a moment-method code that uses the electric field integral 

equation (EFIE) to enforce the boundary condition of the tangential electric 
field being zero at the surface of the antenna of interest. The 
moment-method basis and testing functions used in this code are 
piecewise sinusoidal. 

20 Appendix B lists sample input and output files for the adaptive 

hyperthermia simulation. The first data file was used to generate the 
E-field results for a four auxiliary probe system, and the second data file 
was used to generate the E-field results for a two auxiliary probe system 
The corresponding output files give the values for the array mutual 
25 coupling, quiescent and adaptive transmit weights, channel correlation 
matrix, eigenvalues, and cancellation. 

Referring to FIG. 13, there is shown the hyperthermia phased-array 
antenna system 200 of FIG. 7, redrawn to simplify the following method 
of moments analysis. The RF source 208, power divider 209 and weights 
30 210 n of FIG. 7, are modeled as a plurality of RF signal generators 250, 
through 250 N , feeding its corresponding transmit antenna element 204, 
through 204 N . Each generator 250, through 250 N has a corresponding 
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amplitude and phase weight denoted by w, through w N , and a known 
output impedance Z L . The Ah probe 226 (i.e., the same as the /th probe 
226 of FIG. 7, with different notation) is modeled as a dipole antenna 
having an overall length L p and an open-circuit voltage vf c - induced by 
5 the RF energy transmitted from the antenna array 200. 

The open-circuit voltage at the/th probe antenna 226 is computed 
from the array terminal currents and from Z> the open-circuit mutual 
impedance between the nth array element and the/th probe antenna. Let 
v n°'°' represent the open-circuit voltage at the /th probe due to the nth 
10 transmit-array element. Here, the /th probe can denote either the focal 
point calibration probe (calibration probe 212 of FIG. 7) or one of the 
auxiliary probes used to null a sidelobe. The number of auxiliary probes is 
denoted by N aux . 

Referring also to FIG. 14, the/th probe 226 is modeled as a voltage 
source 260, having an output voltage vf c - , driving a first impedance 
262, representing the input impedance Z /yv of the/th probe, in series with a 
second impedance 264. representing the termination impedance Z r of the 
/th probe. The /th probe receive current tf* flows through these two 
impedances. The output voltage of the/th probe vf ° appears across the 
20 termination impedance Z r 

Referring again to FIG. 13, next, let Z denote the open-circuit 
mutual impedance matrix (with dimensions N X N for the /V-element array). 
The open-circuit mutual impedance between array elements 204 m and 
204 n is denoted Z m>n . It is assumed that multiple interaction between the 
25 hyperthermia array and the auxiliary probe can be neglected. Thus, the 
hyperthermia array terminal current vector / can be computed in terms of 
the transmit weights w as 



l-\Z+Ztfw (3 6 ) 
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Next, let Z n [ be the open-circuit mutual impedance between the Ah probe 
and the nth array element. The induced open-circuit voltage v n ° f c - at the 
fth receive probe, due to the nth array element transmit current can 
then be expressed as 



vtf=ZX (37) 



5 In matrix form, the induced open-circuit probe-voltage matrix v^™- is 



^=*An**-V (38) 



or 



^^b,-Z pf0bmtarra ^Z^Z L tf' 1 w (39) 

where Z pnb9rmny is a rectangular matrix of order N aux X N for the 
open-circuit mutual impedance between the probe array and the 
hyperthermia array. Note that the Ah row of the matrix Z ^ is 

10 written as (ZL zi • . ♦ zJj where i-1 ? */ 

* wnerey- 1,2, • • N aux . The receive 

voltage matrix is then computed by the receiving circuit equivalence 
theorem for an antenna. The receive-antenna equivalent circuit is depicted 
in FIG. 14, where it is readily determined that 



(40) 
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where Z in is the input impedance of the probe. It should be noted that the 
"pro* matrix 's a column vector of length N MX and if c is the fth element 
of the matrix. The probe-receive current matrix is given by 



§nc aft 1 

•orobm- rant.— — (4 -J J 



The fth element of the column vector is denoted if*. /« f,2, • . . 

5 "but 

Finally, the power received by the/th probe is 



where Re means real part. Substituting Equations (40) and {41} into 
Equation (42) yields 



10 



The total interference power received by the auxiliary probe array is given 
by 



P^E/f (44) 
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The incident electric field E is related to the open-circuit voltage 
by the effective height h of the probe antenna as 



,ae, 



hE (45) 



If the length L p of the probe antenna 226 is approximately 0. U or less, 
the current distribution is triangular and the effective height \sh = 0.5L . 
Thus, for a short-dipole probe the open-circuit voltage can be expressed as 



v**Jze {46) 



It then follows from Equation (46) that the £ field for a short-dipole probe 
at position (x,y,z) is given by 



ffcr,2) ^ ft y.*4 (47) 

HP 



finally, the quiescent and adapted E-field radiation patterns are computed 
using the quiescent and adapted weight vectors w q and w„ respectively 
10 In Equations (39, and (47). The moment-method expansion and testing ' 
functions are assumed to be sinusoidal. The open-circuit mutual 
impedances in Equation (39) between thin-wire dipoles in a homogeneous 
conducting medium are computed based on subroutines from the 
moment-method computer code developed by J.H. Richmond In 
1 5 evaluating Z* for the /th auxiliary probe, double precision computations are 
used. 
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As mentioned previously, the array is calibrated (phased focused) 
initially using a short dipole at the focal point. To accomplish this 
numerically, having computed v fo £f , the transmit array weight vector w 
will have its phase commands set equal to the conjugate of the 
5 corresponding phases in v fo * . Transmit antenna radiation patterns are 
obtained by scanning (moving) a dipole probe with half-length / in the 
near-field and computing the receive probe-voltage response. 

The received voltage matrix for the Jth probe (denoted vf) is 
computed at K frequencies across the nulling bandwidth. Thus, 
10 v^ffylv^ff^, . . -.^(frf are needed. For the purposes of this 

computer simulation, the impedance matrix is computed at K frequencies 
and is inverted AT times. The probe channel correlation matrix elements are 
computed by evaluating Equation (2) numerically, using Simpson's rule 
numerical integration. For multiple auxiliary probes, the channel 
T5 correlation matrix is evaluated using Equation (4). Adaptive array radiation 
patterns are computed by superimposing the quiescent radiation pattern 
with the weighted sum of auxiliary-channel-received voltages. 
Wave Pronanation in rnn ^ ct in Q Marih 

To gain insight into the effect of a lossy medium, e.g.. the target 
20 body, on the propagation of an electromagnetic wave, it is useful to 
review certain fundamental equations which govern the field 
characteristics. In a conducting medium. Maxwell's curl equations in 
time-harmonic form are 



VxH=J+/ueE (48) 



and 



Vx£=-/ B ^// (49) 
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where £ and H are the electric and magnetic fields, respectively. J is the 
conduction current density, w=2^is the radian frequency, € is the 
permittivity of the medium, and u is the permeability of the medium. The 
permittivity is expressed as e = ef 0 , where c, is the dielectric constant 
5 (relative permittivity) and c Q is the permittivity of free space. Similarly, p = 
»jj 0 , where p r is the relative permeability andp Q is the permeability of free 
space. For a medium with electrical conductivity a, J and Fare related as 

J*oE (50) 
Substituting Equation (50) into Equation (48) yields 



Vx//=(o+y»c)£ (51) 



From Equations (48) and (49), the vector wave equation in terms of E is 
10 derived as 



&E-y z E=0 (52) 

It is readily shown that 



T =±v>|»(o+/»€) = ± / Wv ^iJ 1 -y_£_ (53) 

N we 



The quantity a/a/e is referred to as the loss tangent. It is common to 
express the complex propagation constant as 
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y-« + yiP (54) 



where a is the attenuation constant and fi is the phase constant. The 
constants a and fi are found by setting Equation (53) equal to Equation 
(54) and then squaring both sides, equating the reai and imaginary parts, 
and solving the pair of simultaneous equations, with the result 



(55) 



5 and 



(56) 



The wavelength A in the lossy dielectric is then computed from 

2n 



(57) 



The intrinsic wave impedance n is given by 




(58) 
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The instantaneous power density of the electromagnetic field is 
given by Poynting's vector, denoted P, 



P -\**H- (59) 



which has units of (W/m*). The time-average power flow density is equal 
to the real part of the complex Poynting's vector. The time-average power 
dissipation per unit volume P a (W/m*) is derived from Maxwell's 
equations, with the result 



v5i-j-4.ui* (60) 



2 



The specific absorption rate (SAR, is the power dissipated or absorbed per 
unit mass (W/kg) of the medium (tissue), or 



161, 



10 



where p is the density of the medium in kg/m* 

It is convenient to have a simple equation for computing the 
propagation loss between any two points in the near field of an isolated 
transmitting antenna. Thus, mutual coupling effects are ignored for the 
time being. Consider a time-harmonic source radiating a spherical wave 
into an infinite homogeneous conducting medium. For an isotropic 
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radiator, and suppressing the e** time dependence, the electric field as a 
function of range r can be expressed as 



£W=£«,-^ (62) 



where E Q is a constant. 

For a source at the origin, the amplitude of the electric field at range 
5 r f is given by 



\BM\-SA^ (63) 



and at range r 2 by 



\^\=E a ±p- ( 64) 



The total propagation loss between ranges r, and r 2 is found by taking the 
ratio of Equations (64) and (63), or 



jjjteHjjg-fe-r.) 



The field attenuation >4 C in dB from range r, to range r 2 due to the lossy 
10 dielectric is simply 
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^-20IOfl 10 («-^>) (66) 



Similarly, the 1/r attenuation loss A r in dB is 



>4,=20log 1( A (67) 

r 2 



Thermal Modeling of an lnhnm 9t)f n« > OU s TW nf T 

A thermal analysis computer program called the transient thermal 

analyzer (TTA), developed by Arthur D. Little, Inc.. has been used to 
5 accomplish the thermal modeling of homogeneous muscle tissue 

surrounded by a constant-temperature water bolus. 

The TTA program uses the finite-difference technique to solve a set 

of nonlinear energy balance equations. Consider a system of 

interconnected nodes that model an inhomogeneous volume for which the 
10 temperature 7} of the Ah node is to be determined. The heat-balance 

equation, which is solved by TTA, Is expressed as 



fo</-^^/5=0 (68) 
M at 



15 



where Q fJ is the net outward heat flow from node / in the direction of 
node/, Pft} is the power into node / at time t. and M, is the thermal mass 
(mass times specific heat) of node /. 

FIG. 15 shows an electric circuit analog 400 which is used to model 
the two-dimensional thermal characteristics of the material volume 402 
which simulates the target body as a plurality of uniformly distributed 
nodes 406 spaced A/ apart. With reference to the Ah node 406,, but 



SUBSTITUTE SHEET 



WO 93/00132 PCT/US92/05464 




m 

applying generally to the other nodes, power P, in watts is delivered 404; 
to the Ah node. Capacitor 408,, having thermal capacitance denoted C f ' 
(with units Joules/°C), is used to model the thermal capacitance at the Ah 
node. Resistor 410 s , having heat resistance denoted R (J (with units 
°C/W), is used to model the heat resistance between Ah node 406,- and 
the/th node 406j. 

With a spacing of LI between nodes (assuming cubic cells), the 
values of R ijt c f and /»,. are computed as 



where k f/ is the thermal conductivity (with units W/m«C) between nodes / 
10 and /; 



Cr?P^ (70) 



where C pl is the specific heat at the /th node and p, is the density (kg/m*) 
at the Ah node; and 



PriSARifiH^ (71) 
where (SARJ, is the SAR for the Ah node, which is given by 

(SAQrgW (72, 



SUBSTITUTE SHEET 



" WO 93/00132 

— PCI7US92/05464 




where a, is the electrical conductivity of the Ah node and |£.| is the 
magnitude of the electric field delivered by the hyperthermia array to the 
Ah node. It should be noted that in substituting Equation (72) into 
Equation (71,, the density P/ cancels. Thus, an equivalent approach to 
computing the power delivered to the /th node is written in terms of the 
t.me-average power dissipated per unit volume of the /th node (denoted 
Paji as 



PrPJtAf (73, 



FIG. 16 is a block diagram showing how TTA is used in the 
,n ^ hermlaSlmu,ationdesc "^ herein. First, the method of moments 
10 500, controlled by the SM. nulling algorithm 502, is used to compute the 
electric field radiation pattern throughout a homogeneous region 
s-mulating muscle tissue, inside an annular phased array 501. These 
E-f,eld simulations assume that the signal received by a short-dipole probe 
wrth.n the region is due to a transmitting phased array embedded in an 
1 5 Infinite homogeneous lossy dielectric (muscle tissue). 

The resulting E-Field power distribution is then read 504 into the 
TTA program 506, which computes the temperature distribution inside an 
e»..pt,cal musc.e-tissue target surrounded with a constant-temperature 
water bolus 507. Because the RF wavelengths in the target and water 
20 bolus are similar, the E-fie»d simu.ations are believed to give a reasonable 
approximation to the field distribution inside the elliptical target The 
computed temperature distribution is output 508 from the TTA for further 
analysis or display. 

The E-fie,d calculation in the assumed infinite homogeneous medium 

Z " T mnal fi6,d attenU8ti0n ° 0t Pr6Sent 10 8 *« h VPerthermia 
system w,th an annular array transmitting through a water bolus into a 

pat,ent. As mentioned earlier, the water bolus has very little RF 

propagation loss. In addition, the transmit array weights are normalized 
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according to Equation (6). Thus, no attempt is made to compute the 
absolute E-field strength in volts/meter in the elliptical target. Instead, the 
peak power in the elliptical target is adjusted (by a scale factor) to produce 
a desired maximum focal-region temperature lT max ) after t minutes. It 
5 should be noted that an approximate absolute scale factor could be 
computed by making an initial computer simulation with an infinite 
homogeneous water bolus and then matching the target boundary field to 
the infinite homogeneous muscle tissue simulation. 

The computer simulation model is related, in part, to the 
tO hyperthermia annular phased-array antenna system shown in FIG. 1 . The 
simulated array is assumed to have a 60-cm array diameter with eight 
uniformly spaced dipole elements which operate over the frequency band 
60-120 MHz. The eight elements of the array are assumed to be fully 
adaptive, whereby seven independent nulls can be formed while 
T5 simultaneously focusing oh a tumor. 

It is further assumed for the purpose of this simulation that the 
adaptive radiation pattern null-width characteristics in a homogeneous 
target are similar to the characteristics observed in an inhomogeneous 
target. The null-width characteristics are directly related to the RF 
20 wavelength, and, only a 5 percent change in wavelength occurs between 
the assumed muscle tissue and water bolus. With this assumption, the 
transmit array may be simulated as embedded in homogeneous tissue, 
which allows direct use of the thin-wire moment-method formulation 
discussed above. 

25 After computing the two-dimensional E-field distribution in the 

homogeneous medium, we then consider only an elliptical portion of the 
homogeneous region and use the ellipse as the homogeneous target. In 
the thermal analysis, the elliptical target is surrounded with a constant 
10°C water bolus. The E-field amplitude is scaled to produce a 46»C 

30 peak temperature, at time t =20 minutes, at the center of the elliptical 
phantom. The initial temperature of the phantom is assumed to be 25 °C 
(room temperature). 
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All computer simulations assume a 120 MHz operating frequency 
with initially four auxiliary nulling probes, i.e., N MX =4. The parameters 
used in the electrical and thermal analyses are summarized in Table 1. 
These parameters are for a frequency of 100 MHz, but is assumed that 
5 similar values of the parameters will exist at 120 MHz. It should be noted 
that although the relative dielectric constants of phantom muscle tissue 
and distilled water are very similar, the electrical conductivities are vastly 
different. The relevant thermal characteristics-density, specific heat, and 
thermal conductivity-are very similar for phantom muscle tissue and 
10 distilled water. 



15 



SIMULATION BESLUJS 

Electric Field for Array i n Hpmnn»n» 91 fff T j ff friff 
Substituting the values /=120 MHz, v=0.5 S/m, and * r =73.5 into 
Equation (53) yields y m = 10.0 + J23.8 for the muscle tissue. With fi m = 
23.8 radians/m, the wavelength in the phantom muscle tissue is A m =26.5 
cm. The attenuation constant for the muscle tissue is a m =10.0 " 
radians/m. Similarly, for distilled water y w = 0.0021 + /22. 5 , so the 
wavelength is A w =27.9 cm. The attenuation constant for the distilled 
water medium is a w =0.0021 radians/m. The propagation loss in the 
20 phantom muscle tissue is 20/og^^ Qf ^ 8? dB/cm ^ 

propagation loss in the distilled water is found to be -0.0002 dB/cm. 
Thus, the total loss due to propagation through 15 cm of distilled water is 
0.003 dB. For 15 cm of muscle tissue the corresponding loss is 13.1 dB 
The wave impedance in the muscle tissue is computed from Equation (58) 
25 as n m = 33.9 + J14.2 Q, and similarly in the distilled water n = 42 1 + 
J0.004 Q. w 

FIG. 17 shows the geometry used in the simulations, which parallels 
the array shown in FIG. 3. A 60-cm-diameter ring phased array applicator 
102 of eight perfectly conducting center-fed dipoles, 104, through 104 8 
30 uniformly surrounds a fictitious elliptical target zone 1 06 with major axis ' 
30 cm and minor axis 20 cm. The length of each dipole array element 
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104 n at 120 MHz in the infinite homogeneous muscle tissue is>U2, or 
13.25 cm. The array focus 107 is assumed at the origin (*=0, yL 0 , 
*=0) and four auxiliary short-dipole probes, 112, through 1 12 4 , with 
length 1.27 cm (0.05,1) are positioned at (x rYf z) coordinates at (15 cm 0 
5 0), (-15 cm, O. 0), (0, 0, 10 cm), and (0, 0, -10 cm), respectively, i.e.,' the 
auxiliary E-field probes are located every 90' in azimuth on the perimeter 
of the target. In rectangular coordinates, each dipole is oriented along the 
9 direction and the feed terminals of each dipole are located at y=0. 
The moment-method computer simulations were run on a Sun 
10 3/260 workstation. The total CPU time for a complete moment-method 
run «s 19.2 minutes. This CPU time includes computing the quiescent and 
adaptive radiation patterns on a 41 by 41 grid of points. The CPU time 
without radiation pattern calculations is 33 seconds. 

FIG. 18 shows the two-dimensional radiation pattern in the plane 
t5 y=0, before nulling, at 120 MHz with uniform amplitude and phase 
illumination. The calculated data are collected on a 41 by 41 grid of 
points over a square region, with side length 76.2 cm, centered at the 
focus 107. The spacing between data points is 1.905 cm, or 0.072A and 
the contour levels are displayed in 10-dB steps. The E-field data are 
20 computed for the case of a 1 .27-cm short-dipole observation probe The 
posmons of the eight dipole radiators 104, through 104 8 are clearly 
evident by the -20 dB contours surrounding each element. The radiation 
pattern is symmetric because of the symmetry of the array and the 
assumed homogeneous medium. 
25 FIG. 19 shows finer contour levels (1-dB steps) for the quiescent 

radiation pattern of FIG. 18. Here, it is evident that the focused main 
beam of the ring array is increasing in amplitude as the observation point 
moves closer to the focus. Away from the main beam region, the pattern 
amplitude is seen to increase as the observation position moves toward 
30 the array perimeter. 

FIG. 20 shows the quiescent radiation pattern of FIG. 18 cut at 
z=0. The large amplitude that occurs at ±30 cm, i.e.. at the position of 
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the phased array applicator 102, is due to the E-field probe's close 
proximity to the transmitting elements 104, and 104 5 . The large 
attenuation that occurs from the array diameter to the focus is due to the 
1/r attenuation loss and the loss in the uniform homogeneous muscle 
5 tissue. FIG. 21 shows the radiation pattern of FIG. 18 cut at x=0. Here, 
the pattern is identical to the pattern of FIG. 20 due to the symmetry of 
the array. In both FIGS. 20 and 21 the boundary of the fictitious elliptical 
target zone 106 is indicated. The target zone of FIG. 20 is larger than that 
of FIG. 21 since the major axis of elliptical target 106 lies along the x-axis, 
10 and the minor axis of target 106 lies along the z-axis. 

The increasing radiation pattern amplitude near the left and right 
sides of the elliptical target of FIG. 20 is shown to produce hot spots in 
the thermal distribution. Because the top (anterior) and bottom (posterior) 
of the elliptical target of FIG. 21 are not as strongly illuminated as to the 
1 5 left and right sides of the elliptical target of FIG. 20, no quiescent hot 
spots occur at the top or bottom. 

Further, FIG. 20 shows that the ring-array half-power beamwidth in 
the target region is approximately 13 cm, or approximately one-half the 
wavelength (26.5 cm) in the phantom muscle tissue. The adaptive nulling 
20 resolution or closest allowed spacing between a strong adaptive null and 
the main beam has been shown to be equal to the half-power beamwidth 
of the antenna. Thus, the closest allowed null position is 13 cm from the 
focus. Since the target width is 30 cm across the major axis, two nulls 
can be formed at <*= ± 15 cm. *=0) at the left and right side of the target 
25 without disturbing the focus. However, if two strong nulls are formed at 
the posterior and anterior (* = 0, z= ± W cm) of the target the focus will 
be compromised. In practice, the water bolus surrounding the target 
would restrict the placement of short-dipole probes 1 12 n to the surface of 
the target. Thus, only weak nulls can be formed at (*-<* *= ± 10 cm) so 
30 that the focus will not be affected by the adaptive nulling process. That 
is. the effect of the two minor axis nulls is to keep the z= ± 1 o cm E-field 
from increasing beyond the quiescent values. 
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Next, adaptive radiation patterns are computed with four auxiliary 
dipole probes 112, through 112 4 positioned as shown in FIG. 17. The 
value of the receiving gain for auxiliary dipole probes 112, and 1 12 2 is 
adjusted to produce a SNR > 35 dB. This amount of SNR results in 
5 greater than 35 dB of nulling in the direction of auxiliary dipole probes 
112, and 1 12 2 . In contrast, the gain values for auxiliary dipole probes 
112 3 and 1 12 4 are turned down to produce about a 3 dB SNR. Thus, 
only about 3 dB of nulling will occur at probe positions 1 12 3 and 1 12 4 as 
the adaptive algorithm reduces the interference to the noise level of the 
10 receiver. The reason for choosing these null strengths will become 
apparent with the data that follow. 

FIG. 22 shows the two-dimensional radiation pattern after nulling 
with four auxiliary probes 112, through 1 12 4 . Two strong adaptive nulls 
at x- ± 15 cm occur as expected, and weak nulling occurs at ± W cm, 
15 also as expected. 

The two strong nulls in the z=0 cut are quantified in FIG. 23, 
where greater than 35 dB of interference nulling or pattern reduction 
occurs at *= ± 15 cm. The peak level at the focus 107 is adjusted to 0 
dB for both the quiescent and adaptive patterns. Two weak adaptive nulls 
20 are seen in the *=0 radiation pattern cut shown in FIG. 24. The weak 
nulls in effect in the adaptive patterns reduce variation from the quiescent 
radiation pattern. 

FIG. 25(a) shows the two-dimensional radiation pattern before and 
after nulling taken longitudinally along a line parallel to the y-axis and 

25 passing through probe 112, as shown in FIG. 25{b), i.e., x=15 cm, z=0. 
This radiation pattern clearly shows that a strong adaptive null also 
extends in the y direction from the E-field probe being nulled. 

FIG. 26(a) shows the two-dimensional radiation pattern before and 
after nulling taken along the y-axis and passing through the focus 107 as 

30 shown in FIG. 26(b). i.e., x=0, 2 -0. This radiation pattern clearly shows 
that the E-field at the focus remains virtually the same in the y direction 
before and after adaptive nulling at the E-field probes 112, through 1 12 4 . 
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RG. 27(a) shows the transmit array amplitude weights before (solid 
line) and after (broken line) nulling, and FIG. 27(b) shows the transmit 
array phase weights before (solid line) and after (broken line) nulling. As 
shown, the adaptive transmit weights exhibit a 5-dB dynamic range in FIG 
5 27(a). 

FIG. 28 shows the channel correlation matrix eigenvalues before 
(solid line) and after (broken line) nulling. There are two large eigenvalues, 
A, and A 2 , and two weak (non-zero) eigenvalues, A z and A 4 , shown in FIG. 
28. These eigenvalues are directly associated with the two high-SNR 

10 auxiliary probes 112, and 1 12 2 , and the two weak-SNR auxiliary probes 
1 12 3 and 1 12 4 , respectively. Note that the 0-dB level in FIG. 28 is equal 
to the receiver noise level. The probe-array output power before and after 
adaptive nulling is 31.4 dB and 0.9 dB, respectively, as calculated from 
equation (10). This difference in power before and after nulling indicates 

15 that the adaptive cancellation is -30.5 dB. 

Temperature r)istrih.,tinn In Ff liPtical Ph*n fm 

To simulate the temperature distribution in the target body resulting 
from the calculated E-fields. the transient thermal analysis (TTA) software 
is used to compute the temperature distribution in an elliptical phantom 
20 surrounded with a constant-temperature water bolus. The 41 X 41 

two-dimensional E-field radiation pattern data of FIGS. 18 through 24 are 
used as the power source for the thermal node network. Two node 
spacings are considered. First, the node spacing Ax . Az «= A/ = 1.905 
cm (coarse grid) is used to obtain thermal data. Then, the node spacing is 
25 decreased by a factor of two to LI = 0.9525 cm (fine grid) to check 
convergence. The coarser spacing is shown to be adequate. 

The scale factors used to convert the normalized E-field distributions 
to a power level that induces a 46*C peak temperature at t =20 minutes 
are 94.1 dB and 96.0 dB for the quiescent and adaptive patterns, 
30 respectively. These scale factors are determined empirically. From 

Equations (68) through (73) and the parameter values given in Table 1, all 
resistors R fJ in the phantom muscle tissue had a value of 96.5«C/W and 
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all resistors R u in the water bolus had a value of 87.2°C/W. The value of 
the capacitors C t In the phantom muscle tissue is 23.6 J/°C. Capacitors 
are not used in the water-bolus region in the input to the transient thermal 
analysis software. Instead, a constant temperature of 10°C is enforced at 
5 each water-bolus node. With a 41 X 41 grid, a total of 3280 resistors and 
1681 capacitors are used in the thermal simulation. The CPU time 
required to compute this temperature distribution is under four minutes. 

FIG. 29 shows the two-dimensional temperature distribution 
produced at time t =20 minutes in the elliptical phantom muscle tissue 
10 target 106 without adaptive nulling. To generate FIG. 29, the power 
source used in the transient thermal analysis is the quiescent radiation 
pattern given in FIG. 18. The initial temperature (at time t=0) is 25°C. 
Notice the occurrence of two hot spots 122 and 124 on the left and right 
sides of the elliptical phantom, respectively. The peak temperature at 
15 focus 107 is 46«C, which is achieved by scaling the normalized quiescent 
E-field as described earlier. The two hot spots 122 and 124 are quantified 
in the z=0 temperature pattern cut shown in FIG. 30, and have a peak 
temperature at each hot spot of approximately 41 *C. The temperature 
profile for x=0 in FIG. 31 shows no hot spots. As any undesired hot spot 
20 is a potential source for compromising the therapy session, adaptive 
nulling is used to reduce the sidelobes corresponding to the hot spots. 

FIG. 32 shows the simulated two-dimensional thermal distribution at 
time t =20 minutes, with adaptive nulling at four auxiliary probes 112, 
through 1 12 4 in effect. The focal-spot diameter at focus 107 with 
25 adaptive nulling is equivalent to the focal-spot diameter before adaptive 
nulling, shown in FIG. 29. Hot spots on the left and right sides of the 
target 106 are eliminated. FIG. 33 shows a comparison of the 
temperature distribution before (solid line) and after (broken line) nulling 
along the major axis (z=0) of the target ellipse 106. Similarly, FIG. 34 
30 shows the temperature distribution before (solid line) and after (broken 
line) nulling along the minor axis (x=0) of the target ellipse 106. 
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The convergence of the previous thermel simulations was verified 
by increasing the density of E-field observation probe positions by a factor 
of two, with a new spacing between points of 0.9525 cm, still with a 41 
X 41 grid. The ring array operates as before at 120 MHz, and there are 
5 four auxiliary probes 112, through 112 4 laid out as shown in FIG. 17. As 
the auxiliary positions are the same, the adaptive weights and channel 
correlation matrix eigenvalues in FIGS. 27 and 28, respectively, remain the 
same. From the parameter values in Table 1. all resistors R tj in the 
finer-grid muscle-tissue phantom had a value of 193.0«C/W and all 
TO resistors R tJ in the water bolus had a value of 1 74.4-C/W. The value of 
the capacitors C, in the phantom muscle-tissue is 2.95 J/°C. Again, a 
constant temperature of 10°C is enforced at each water-bolus node'. The 
E-field scaling factors to raise the focal-point temperature to 46«C before 
and after nul.ing are 76.5 dB and 78.4 dB, respectively. The finer-grid 
1 5 two-dimensional thermal distributions before and after nulling are shown in 
FIGS. 35 and 34. respectively. Although the temperature contours are 
smoother, the general agreement between these patterns and the 
coarser-grid patterns in FIGS. 29 and 32 are evident. Similarly, 
one-dimensional thermal pattern cuts with the finer grid are shown in 
20 FIGS. 37 <* axis) and 38 U axis), and good agreement with the coarse-grid 
patterns of FIGS. 33 and 35, respectively, is observed. In particular, the 
finer detail in FIG. 37 shows that the hot spots 122 and 124 are at 
approximately 42"C compared to 41 « C observed for the coarse grid of 
FIG. 33. Thus, convergence of the coarse-grid thermal patterns is 
25 demonstrated. 

Elliptical Arrfly 

An elliptical phased-array hyperthermia applicator, having a 70 cm 
major axis and a 60 cm minor axis, was also analyzed by computer 
simulation. The computer simulation parameters were the same as those 
30 applied to the analysis of the annular array. Generally, the computer 

simulations show that reduced hot spot temperatures are observed along 
the major axis of the elliptical phantom, without adaptive nulling, while 



SUBSTITUTE SHEET 



WO 93/00132 PCT/US92/05464 

* m 

small increases in hot spot temperatures occur along the minor axis. 
Certain tumor geometries may be heated more efficiently with an elliptical 
array than with an annular array. 
EXPERIMFMTAt BfSULIS 
5 Experimental data have been gathered from a commercial annular 

phased-array hyperthermia system modified to perform a gradient search 
algorithm to produce an adaptive null (or focus) at one or more auxiliary E- 
field probe positions. The results confirm that a strong null can be formed 
at the surface of the target body without significantly affecting the power 
1 0 delivered at the focus of the hyperthermia system. 

FIG. 39 shows a cross-sectional view of the experimental system 
700, which is a modified BSD-2000 SIGMA-60 annular phased-array 
hyperthermia applicator, available from BSD Medical Corporation, Salt Lake 
City, Utah. The annular array antenna 702 of the system is 59 cm in 
1 5 diameter and includes eight uniformly spaced dipole antennas 704, 

through 704 8 , excited with a four channel transmitter at 100 MHz. Each 
of the four transmit channel signals are distributed by separate coaxial 
cables from the hyperthermia controller (not shown) to a two-way power 
divider having two outputs. The two outputs of each two-way power 
20 divider drive a pair of dipole antenna elements through a pair of coaxial 
cables. 

A cylindrical phantom target body 706 is supported by a patient 
sling 705 which centrally locates the phantom within the annular array so 
that the longitudinal axes of the phantom and the annular array correspond 
25 to each other. Phantom target 706 is a 28 cm diameter X 40 cm long 
polymer bottle filled with saline solution (0.9% NaCI), which simulates a 
human subject. A deionized water-filled bladder 705 provides a water 
bolus between the annular array and the target phantom. 
Sinole Artanflx/p m.,ii 

30 Three E-field probes are used to monitor the amplitude of the E-field 

at various sites in and around the phantom for this experiment. The first 
E-field probe 715 (BSD Medical Corp. Model EP-500) is located inside the 
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phantom at the center, or focus. 707 of the array which simulates the 
tumor site. This probe monitors the amplitude of the E-field at the tumor 
site as the null is formed at the null site. The second E-field probe 720 
(BSD Medical Corp. Model EP-100) is taped onto the outside surface of 
5 the phantom at the desired null location which simulates an E-field probe 
taped to a patient's skin. The probe is used to monitor the amplitude of 
the E-field at the null site as the null is formed by the gradient search 
algorithm. The third E-field probe 721 (BSD Medical Corp. Model EP-400) 
is taped onto the outside surface of the phantom diametrically opposite 
10 the location of the null site probe 720. This probe is used to monitor the 
amplitude of the E-field away from the null and focus sites and provides an 
E-field amplitude reference for the experiment. 

The transmit array amplitude and phase control software and the 
electric field probe monitoring software supplied with the BSD-2000 
1 5 system were modified to incorporate a gradient search feedback routine 
for adaptive nulling and adaptive focusing. Pascal source code listings and 
sample output of the adaptive nulling and focusing gradient search 
feedback routines are attached hereto as Appendicies C and D, 
respectively. 

Fig. 40 graphically illustrates the results of this experiment, showing 
the measured E-field probe amplitude, in dB. versus the gradient search 
iteration number. The dB values are obtained by computing 10log w (probe 
output signal) and normalizing the resulting values to O dB at iteration 0. 
It is evident from this graph that the gradient search formed a strong E- 
25 field amplitude null at the null site, on the order of -1 5 to -20 dB with 

respect to the reference site, in less than 50 iterations. (The apparent rise 
in the E-field amplitude at the null site between iterations 45 and 50 is 
most likely due to noise associated with the convergence calculations). 
Furthermore, the measured E-field amplitude at the tumor site was reduced 
by no more than -5 dB with respect to the initial reference level. Fig. 41 
shows a graphic illustration of the power, in dB. calculated at the null site 
versus gradient search iteration. It is evident from this graph that the 
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gradient search causes the null site power to monotonlcally decrease with 
each iteration, achieving an approximately 12 dB reduction in power within 
50 iterations. 

It should be noted that for at (east the first 30 iterations of the 
5 gradient search, there is good agreement between the computer 
simulations, presented above, and these experimental measurements. 
After approximately 30 iterations, however, the results of the computer 
simulations differ from the experimental measurements. One reason for 
this difference is that the computer simulations herein described do not 
TO attempt to accurately model all the characteristics of the BSD-2000 
system used for the experiments. For example, the simulations do not 
account for phase shifter non-linearities. A/D convenor errors, or D/A 
converter errors associated with the system which will affect the 
experimental measurements, especially at the relatively low signal levels 
15 present after 30 iterations. Thus, it is not expected that the computer 
simulations and the experimental results will necessarily behave the same 
where the signals or computations are most affected by measurement 
system noise. 

Single Adantive N.iH. iw ph p „ t — 
20 Referring to Fig. 42, in another experiment, a single adaptive null 

was produced in a beef phantom 706' used in place of the saline phantom 
706 of Fig. 39 to better simulate human tissue. Beef phantom 706' was a 
24 lbs. hind leg cut having a front face width of 38 cm, a front face 
height of 23 cm and a thickness of 15 cm. The E-field probe positions 
25 used with the beef phantom are analogous to the E-field probe positions 
used with the saline phantom. That is, the first E-field probe 715' (BSD 
Medical Corp. Model EP-500) is located inside the beef phantom at the 
center, or focus, 707' of the array which simulates the tumor site. This 
probe monitors the amplitude of the E-field at the tumor site as the null is 
30 formed at the null site. The second E-field probe 720' (BSD Medical Corp. 
Model EP-100) is taped onto the outside surface of the beef phantom at 
the desired null location which simulates an E-field probe taped to a 
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patient's skin. This probe is used to monitor the amplitude of the E-field 
at the null site as the null is formed by the gradient search algorithm. The 
third E-field probe 72V (BSD Medical Corp. Model EP-400) is taped onto 
the outside surface of the beef phantom diametrically opposite the location 
5 of the null site probe 720'. This probe is used to monitor the amplitude of 
the E-f,eld away from the null and focus sites and provides an E-field 
amplitude reference for the experiment. 

Rg. 43 graphically illustrates the results of the beef phantom 
experiment, showing the measured E-field probe amplitude, in dB, versus 
the gradient search iteration number. Again, the dB values are obtained by 
computing 1 0log 10 iprobe output signal) and normalizing the resulting 
values to 0 dB at iteration 0. It is evident from this graph that the 
gradient search formed a strong E-field amplitude null at the null site 720' 
on the order of -18 to -20 dB with respect to the reference site, in less 
1 5 than 50 iterations. Furthermore, the E-field amplitude at the tumor site 
was reduced by no more than -2 dB with respect to the initial reference 
•evef. It should be noted that these results are very similar to the results 
obtained with the saline phantom (Fig. 40). 

Fig. 44 shows a comparison between the temperature rise at the 
20 beef phantom tumor site 707' and the null site 720' during nulling. A 
thermocouple probe was located at each of the tumor and null sites, and 
the RF power was applied in four intervals of 15 minutes power on and 5 
minutes power off for a total experiment time of 80 minutes. The gradient 
search performed 10 iterations during the 1 5 minute power on portion of 
25 each interval. Temperature measurements were taken during the 5 minute 
power off portion of each interval, i.e.. one measurement for each 10 
Iterates. The experimental data shows that the tumor site was initially 
at about 26°C and the null site was initially at about 27"C before applying 
RF power. After 40 minutes (30 minutes power on and 10 minutes power 
30 off, the temperature of the tumor site has risen 4°C to about 30°C while 
the null site has risen only fc to about 28°C. After 80 minutes (60 
mmutes power on and 20 minutes power off, the temperature of the 
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tumor site has risen 8°C to about 34°C, while the temperature of the null 
site has risen only 3°C to about 30°C. Thus, an approximate differential 
of about 4°C is attained between the tissue temperature of a deep-seated 
target and the temperature of a single surface null site in a beef phantom. 

5 Two AdantivA fl if » ff 

In another experiment, two adaptive nulls on the surface of the 
saline-filled cylindrical phantom were generated and measured. This 
experiment used the same configuration as shown in Fig. 39 with the one 
E-field probe 715 (BSD Medical Corp. Model EP-500) located inside the 
phantom 706 at the center, or focus, 707 of the array simulating the 
tumor site. Two E-field probes 720 and 721 (BSD Medical Corp. Model 
EP-100) were located on the outside surface of the phantom at 
diametrically opposite positions representing the two non-Invasive null 
sites. Probe 715 monitored the E-field amplitude at the tumor site while 
15 probes 720 and 721 monitored the E-field amplitude at the null sites. 

FIG. 45 graphically illustrates the results of this experiment, 
showing the measured E-field probe amplitude, in dB, versus the gradient 
search iteration number. Again, the dB values are obtained by computing 
10log 10 (probe output signal) and normalizing the resulting values to 0 dB 
20 at iteration 0. It is evident from this graph that the gradient search formed 
two strong adaptive E-field amplitude nulls at the null sites, on the order of 
-10 to -20 dB with respect to the reference site, in about 50 iterations in 
particular, at Iteration number 50 the null strength at probe 720 (probe site 
2) is approximately -18.0 dB and the null strength at probe 721 (probe site 
25 3) is approximately -1 1 .8 dB. Furthermore, the E-field amplitude at the 
tumor site 715 (probe site 1) was held close to a constant value (0 dB) 
throughout the 50 iterations. 
Adaptive Foriiiyinq 

In another experiment, adaptive phase focusing was used to 
30 maximize the E-field amplitude at a selected location different from focus 
707 of the saline-filled cylindrical phantom 706 of FIG. 39. In this case 
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the selected focus site was at E-fleld probe 720 (BSD Medical Corp. Model 
EP-100) located on the outside surface of the cylindrical phantom. 

FIG. 46 graphically illustrates the results of this experiment, 
showing the measured E-field probe amplitude, in dB, versus the gradient 
search iteration number for 30 iterations. Again, the dB values are 
obtained by computing fOfog f0 f prob e output signal) and normalizing the 
resulting values to 0 dB at iteration 0. The initial phase weights applied to 
the transmit elements of the array were equal, nominally producing an E- 
field focused at the center of the array 707. The gradient search was 
used to adjust the phases of the array transmit weights to maximize the E- 
field amplitude at probe site 720, i.e., refocus the array at probe 720 The 
transmit weight amplitudes were held constant over the 30 iterations As 
shown, the gradient search converged in about 10 iterations and the 
power at probe 720 increased by about 0.9 dB compared to its initial 
value. This result demonstrates that adaptive focusing can be used 
successfully to optimize the peak power delivered to a tumor site 



Clinical Applicator, 

A modified BSD-2000 Sigma 60 system can be used as a clinical 
adapt.ve hyperthermia system for implementing the adaptive nulling and 
20 focusing techniques of this invention. An unmodified BSD-2000 

hyperthermia system uses four transmit channels to energize the eight 
transmit elements (in pairs) of the annular array, and eight EP-400 (or EP- 
100) non-invasive E-field probes to monitor clinical hyperthermia 
treatments. The eight E-field probes can provide feedback signals to the 
25 controller performing the adaptive nulling and/or focusing algorithms 
Theoretically, three independent adaptive nulls (and/or peaks) can be 
formed by adaptively adjusting the phases and gains of the four transmit 
channels. Any three of the eight E-field probes can provide the feedback 
s,gna«s required to produce a null (or peak, at the corresponding probe 

Various treatment protocols are possible for selecting desired null 
sites, depending on the particular patient and case history. One protocol 
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would place the eight E-field probes around the circumference of the 
patient to measure the E-field strength at each probe before nulling and 
thereby identify the strongest electric fields on the surface of the patient 
indicating potentially serious hot spots. Adaptive nulling would then be 
5 applied to minimize the electric field at the three probes having the 
strongest electric fields before nulling. Alternatively, if the patient can 
localize a painful hot spot during treatment, adaptive nulling would be 
applied to minimize the electric field at the E-field probe closest to the 
identified hot spot. 

10 The number of adaptive nulls required will vary with patient and 

pathology. In some situations it is possible that more than three 
independent adaptive nulls will be required to achieve a therapeutic 
thermal distribution in the patient. In such a case, an extension of the 
four transmit channel BSD-2000 system to an eight transmit channel 
15 configuration will allow up to seven independent adaptive nulls and an 
adaptive focus to be formed. 

MQNOPOI F ARRAY FMRnnm/ifM-r 

FIGS. 47-49 show an embodiment of a non-invasive RF monopoie 
20 phased-array hyperthermia system 400 for treating malignant brain 
tumors. Hyperthermia system 400 features a monopoie phased-array 
transmit antenna 402 having a plurality of monopoie transmit antenna 
elements 404 piaced in proximity to the cranium 406 of a patient to be 
treated for a malignant brain tumor 407. Hyperthermia system 400 also 
25 features improved focusing characteristics through using one or more non- 
invasive electric field probes 412 placed on or near the patient's cranium 
-n conjunction with the near-field adaptive focusing and nulling apparatus 
and methods of this invention. 

^ On0po,e Phased-array 402 is used to therapeutical heat the brain 
30 tumor 407. typical.y located 1-3 cm below the sxin surface of the cranium 
406. by adaptively focusing the RF electric field energy radiated by the 
monopoie radiator elements 404 into the tumor 407. In practice it is 
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dangerous, and often impossible, to invasively place an E-field probe into 
the brain tumor site to facilitate adaptive focusing of the RF energy into 
the tumor. "Hot spots" are not typically a problem with the monopole 
hyperthermia array described since the side lobes generated by a phased- 
5 array in this near-field geometry are much lower than those generated by 
the annular phased-array described above. Thus, the monopole phased- 
array receives its major benefit by applying the adaptive focusing of this 
invention to more precisely focus energy into the tumor site. However, 
the adaptive nulling of this invention may also be used if -hot spots" do 
10 develop through the use of the monopole array. 

The monopole phased array antenna 402 is mounted inside an 
water-tight enclosure 430 having a generally circular top and bottom 
surface 432 and 434. respectively, and a cylindrical or conical side surface 
436 connecting the top and bottom surfaces. The enclosure is made from 
15 non-conductive plastic material, such as plexiglas, but may also be 

fashioned from any material which acts as an electrical insulator and will 
not interfere with the RF radiation patterns generated by the monopole 
phased-array inside the enclosure. 

The bottom surface 434 has a central elliptical aperture 438 which 
20 accommodates a portion of the patient's cranium 406 to allow tumor 407 
to be disposed within the interior of enclosure 430 adjacent to the 
monopole phased-array 402. A flexible silicone rubber membrane 440 
covers the aperture to maintain the water-tight integrity of the enclosure 
The enclosure 430 can be filled with chilled de-ionized water 442 for 
25 cooling the patient's skin during hyperthermia treatment. The de-ionized 
water can be temperature controlled and circulated through the enclosure 
430 to maximize the cooling effect. 

Bottom surface 434 also includes an RF conducting ground plane 
444 mounted co-planar with the bottom surface and which acts as an RF 
30 reflector for monopole antenna elements 404. This ground plane may be 
fashioned out of a metal sheet, metal foil, metal mesh, or any other RF 
conductive materia, which can be fashioned to cover the area of bottom 
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surface 434. The ground plane may also be Imbedded into the bottom 
surface by, for Instance, laminating the ground plane between two layers 
of insulating material. 

The monopole radiator elements 404 are each mounted 
5 perpendicularly onto bottom surface 434 so that they may be energized 
from outside the enclosure, yet remain insulated from the ground plane 
444. In one preferred embodiment, each monopole radiator element 404 
is a 1/4 wavelength long straight wire radiator threadably attached on one 
end into a connector mounted onto surface 434 and insulated from ground 
10 plane 444. It is understood that another form of a monopole antenna 
element, other than a straight wire radiator, can also be used. For 
instance, helical monopole, conical monopole, and sleeve monopole 
antenna elements are also appropriate for use as monopole array elements 
of the present invention. 

15 Each monopole element is energized through a coaxial cable fed 

through the bottom surface 434 to the connector. The bottom surface 
may also be provided with extra monopole connectors which allow 
repositioning of the monopole radiator elements within the enclosure 
Repositioning allows the user to change the geometry of the monopole 
20 phased-array antenna as well as position the antenna adjacent to the 
tumor location to maximize the therapeutic effect. 

An RF reflecting screen 446 (FIGS. 48 and 49} can be placed behind 
the monopole antenna elements 404 to direct more of the radiated RF 
energy toward the cranium, i.e., energy which would otherwise be lost 
through the side of the water-tight enclosure. Reflecting screen 446 is 
typically positioned in the water bolus approximately 1/4 wavelength 
behind the monopole antenna elements and has a cylindrical reflecting 
surface extending perpendicularly from the ground screen to a height of 
approximately twice the wavelength of the radiated energy. Alternatively 
the reflecting surface of the screen can be curved toward the target to 
further enhance radiation of the target. The reflecting screen is 
constructed from high frequency RF conducting mesh which allows water 
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to freely flow through it, and is electrically connected to the ground screen 
using good high frequency RF construction practices. 

Each monopole radiator element 404 is configured as a 1/4 
wavelength radiator to resonate at approximately 915 MHz which is 
5 effective for heating tumors 1 to 3 cm. or more, beneath the surface of 
the patient's skull. The monopole phased-array can include a varying 
number of radiator elements spatially arranged in a variety of patterns. 
The spacing between the monopole antenna elements is typically between 
1/2 to 1 wavelength. Furthermore, the number and location of the non- 
10 invasive electric field probes 412 can also vary depending on the 
hyperthermia focusing patterns desired. 

In the adaptive hyperthermia monopole phased-array of this 
invention, non-invasive E-field probes are used in conjunction with the 
adaptive focusing apparatus and techniques of this invention to maximize 
1 5 the RF power delivered to the tumor site inside the cranium. Computer 
simulations, presented herein, show that the optimum focused (e.g., with 
an invasive E-field probe) phased array can produce an RF energy pattern 
with maximum electric field strength at the tumor site and no undesired 
hot spots within the cranial target. 
20 With the adaptive hyperthermia monopole phased-array described 

herein, RF energy peaks are adaptively formed to maximize the electric 
field energy delivered to the target focus. As shown, the focused energy 
peak achieved by the adaptive focusing apparatus of this invention is 
invasive to the cranial target extending into the tumor region. 
25 Referring to FIG. 50. the adaptive-focusing monopole phased-array 

hyperthermia system of this invention can be described in terms 
corresponding to the generalized annular phased-array system schematic 
diagram of FIG. 6. Specifically, monopole transmit elements 404 n of 
hyperthermia transmitting antenna array 402 correspond respectively to 
30 the dipole transmit elements 104 n of annular phased array applicator 102 
of FIG. 6. Furthermore, the plurality of E-field auxiliary probes 41 2 m . 
correspond to the E-field probes 1 12 m of FIG. 6. It is apparent that'he 
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monopole phased-array hyperthermia system herein described does not 
take advantage of an electric field probe placed at the tumor 407 
analogous to the receiving probe 115 used with the annular phased array 
applicator 102 of FIG. 6, to maximize the focus radiated energy into the 
5 tumor. Use of a probe at the focus wouid in most cases require a surgical 
procedure to invasively place the probe within the patient's brain. 

The receiver 1 14, signal processor 1 16, RF source 108, and 
weighting functions 1 10 n (FIG. 6, operate with the monopole array as 
described above with regard to the annular array, except that the signal 
10 processor 1 1 6 performs an adaptive focusing algorithm described below 
which is related to the adaptive nulling algorithm. That is, signal 
processor 116 performs either a sample matrix inversion (SMI, algorithm 
or a gradient search algorithm on the signals output from receiver 1 14 and 
updates the adaptive array weights w n (with gain g and phase 0, to 
1 5 rapidly (within seconds) focus energy at the tumor 407. 

Referring to HQS. 51 and 52, there is shown an analytical mode, of 
an embodiment of an eight-element, 915 MHz hyperthermia monopole 
Phased-array 402 of FIG. 47. Phased-array 402 has transmit antennas 
404, through 404 8 , arranged adjacent to an elliptical phantom target 406 
20 representing the cross section of the human cranium at the tumor level 
The focus 407 of the elliptica, phantom models the location, approximate* 
2.0 cm below the surface of the cranium, of the brain tumor to receive 
hyperthermia treatment, i.e.. the focus of RF energy for the phased array 

25 tr^dTa r h 0,US ^ " " * ^ ^ «*< - * 
25 treated as a homogeneous medium for analysis purposes 

The monopole radiator elements 404, through 404, are arranged as 
a 1 20 crcular arc array of uniformly spaced elements having a constant 
rad,us of 12.7 cm relative to the geometric center of the cranium C. i e 
at x - 0.0 cm, z = 0.0 cm. The tumor site, or focus 407 of the RF 
30 energy, ls assumed to be at x = 0.0 cm, z - 5 .08 cm for simulation 
purposes. (In an alternative preferred embodiment the monopole array 
elements form a circular arc having a geometric center at focus 407 
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(target) rather than at the center of the cranium C. This has the 
advantage that less phase focusing should be required to maximize the 
energy delivered to the focus, and thus the required number of gradient 
search iterations is reduced.) 
5 Six auxiliary RF E-field probes, or sensors, 412, through 41 2 6 (i.e., 

receiving antennas) are placed on and near the perimeter of the elliptical 
target to model non-invasive E-field probes placed on and near the skin of 
the cranial target. Auxiliary probes 412,, 41 2 2 , and 41 2 3 are uniformly 
spaced in an arc row. between the arc array and the cranial target, having 
10 a constant radius of approximately 4.0 cm relative to the desired focus 
407. Auxiliary probes 41 2 4 , 41 2 5 and 41 2 6 are placed on the target skin 
adjacent to focus 407. Specifically, the first arc row of electric field 
probes 412, through 41 2 3 may be denoted as probes L v M v and R u 
respectively, and the row of electric field probes 41 2 4 through 41 2 6 may 
1 5 be denoted as probes L 2 , M 2 , and R 2 , respectively. The electric field 
probes are arranged so that corresponding probes on the two rows are 
located along a radial line extending from the desired focus at tumor site 
407 and are spaced 1/4 to 1/2 wavelength apart. That is, probe pair 
iL v L 2 ) is located along radial r L . probe pair W,.M 2 ) is located along radial 
20 r M and, probe pair (R V R 2 ) is located along radial r R . 

The gains and phases of the monopole elements are adaptively 
adjusted as described beiow to focus the energy output from the 
monopole phased-array into the tumor site 407 located several centimeters 
below the surface of the cranium. From the phased-array geometry of Fig. 
25 52 it is observed that an electric field focused at tumor 407 will be 

balanced and symmetric with respect to the line x=0. Furthermore, the 
electric field is attenuated in the water bolus external to the cranium in the 
direction away from the transmit array. To achieve a focus interior to the 
cranium at the tumor site 407 it is assumed that the amplitude difference 
30 between the electric field adjacent to the skin surface of the cranium and 
the field approximately one quarter wavelength exterior to the skin surface 
of the cranium must be constrained to a desired value. This desired value 
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»s typically a minimum to avoid "hot spots" on the skin surface. Similarly 
the amplitude of the electric field in the transverse direction should be 
bafanced to minimize the electric field variation between the left and right 
electric field probes with respect to the middle electric field probes I e 
maintain electric field symmetry with respect to the x-axis. The electric 
field differences in the radial direction may be denoted by 

&A L1e =\A u -A u \ , ( 74 j 



AA UT2-\Aui-Ausl . (75) 

LA B1£-\ A fir A f{z\ • (76) 

where A denotes the amplitude of the electric field measured by the 
specified field probe, p is a unit vector in the radial direction which bisects 
the transmit monopoie array and^r is a transverse unit vector as shown in 
tO FIG. 52. The electric field differences in the transverse direction for the 
first row may be denoted by 

*Aun-\Au-Aw\ . (77) 

and, 

AA mn B \Am-Autl . (78) 

and the electric field differences in the transverse direction for the second 
row may be denoted by 
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* a ub b \Au-Avb\ . (79) 

and, 

& A fui2~ lAfv-Aj^l (80) 

A figure of merit F can be defined as 

F=«(*A L12 +AA m2 +AA M2 ) 

+ AA uii*AA mn +AA tMg +AA fii e , (81) 



10 
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where a is a scale factor used to adjust the effect of the eiectrical field 
gradient caused by attenuation in the radial direction between the phased- 
array antenna and the target. The figure of merit F involves seven 
constraints which are easily taken into account by the eight transmit 
element phased-array described. The gradient search algorithm described 
above is used to minimize the figure of merit F. 

FIG. 53 shows a simulated two-dimensional quiescent radiation 
pattern in the plane y=0 lor the eight element monopole arc array of FIG 
52 operating at 915 MHz before adaptive focusing, i.e., with uniform 
amplitude and phase illumination. This radiation pattern was calculated 
using the moment-method described above, and the calculations assume 
an infinite homogeneous conducting medium simulating phantom brain 
tissue, i.e., e r = 50.0, <r=1.3. The focus of the array is atx=0.0, z=5 0 
cm. The positions of the eight monopole radiators 404, through 404 a are 
clearly evident by the -10 dB contours surrounding each element. The 
radiation pattern is symmetric because of the symmetry of the array and 
the assumed homogeneous medium. 

FIG. 54 shows the quiescent radiation pattern of FIG. 53 cut at 
z=5.0 cm, which is through the tumor site 407 assumed to be at *=0.0, 
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*=5.0 cm. The focused main beam is centered at *=0.0 as desired. The 
half-power beamwidth is approximately 2.0 cm, which is close to 1/2 
wavelength for the full ring array. 

FIG. 55 shows a simulated two-dimensional thermal pattern 
5 expected for the quiescent radiation pattern of FIG. 53 at time t =20 
minutes. This simulation assumes that the elliptical phantom brain tissue 
•s surrounded by a 10"C constant temperature water bolus and that the 
initial temperature of the brain tissue phantom is 25 °C. 

FIG. 56 shows a configuration of the monopole phased-array 
10 hyperthermia system of this invention configured to uniformly heat a large 
•ntra-cranial tumor target 407'. In this case, a set of auxiliary E-field 
probes 412 are uniformly spaced along the skin surface of the cranium 
between the monopole radiator elements 404 and the tumor site 407' 
Here, seven E-field probes 412, through 41 2 7 are used, denoted P, 
1 5 through P 7 respectively, but the quantity of probes required will vary 
according to the tumor size and location. The gains and phases of the 
monopole elements are adaotively adjusted to uniformly distribute the 
electric field energy at the E-field probes 412, through 412 7 . From the 
geometry of FIG. 56 it is observed that a uniform electric field at the E- 
20 field probes will produce a substantially uniform electric field inside tumor 
407' to induce uniform heating of the tumor. 

To achieve a uniform electric field distribution across all the E-field 
probes, and thus uniform heating of the tumor 407'. the amplitude 
difference between any two adjacent E-field probes on the skin surface of 
25 the cranium must be minimized. The electric field differences between 
adjacent E-field probes may be written as 

AA Pi2=\A„-A K \ , <82 , 



AA P23-\Ap S -A ps \ , (83) 
AA P34=\Ap 3 -A P4 \ , (84) 
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AA P45=\Ap4-Aps\ . (85) 

A */w=M«-*»l . (86) 

and, 

A VIV^I . (87) 

where >4 denotes the amplitude of the electric field measured by the 
specified field probe. A figure of merit F may be defined as 

F=AA ^ AA ^A^AA^AA^AA^ . (88) 

The figure of merit F involves six constraints on measured differences and 
an additional constraint on total power radiated by the transmit array, all 
of which are easily taken into account by the eight transmit element 
phased-array described. The gradient search algorithm described above is 
used to minimize the figure of merit /^and thereby achieve a uniform 
electric field distribution across the E-field probes. 

It should be noted that this approach for achieving a uniform electric 
f.eld distribution to effect uniform heating within a large mass is not 
bmited to the monopole array heating of a brain tumor, and for example 
can be similarly implemented with the annular phased-array hyperthermia 
applicator of FIG. 1 for uniformly heating a large mass «n other areas of the 
15 body. 

FIGS. 57 and 58A show an alternative preferred embodiment of the 
monopole phased array applicator 400 of FIG. 48 including the addition of 
a top ground plane surface 450 positioned above the monopole antenna 
elements 404 n and extending from the reflecting screen 446 toward the 
20 target body parallel to ground plane surface 442. The top ground plane 
surface 450 combines with ground plane surface 442 to form a parallel 
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plate waveguide region 452 between the monopole antenna elements and 
the target body 406. The spacing between the parallel plates {i.e., 
between surfaces 442 and 450) can be used to adjust the radiation 
pattern in the direction perpendicular to the parallel plates. The spacing 
5 between the parallel plates is typically between 1/2 and 5 wavelengths. 
FIGS. 58B and 58C show alternative preferred embodiments of the parallel 
plate waveguide of FIG. 58A having non-parallel waveguide surfaces, and 
flared waveguide surface forming a horn, respectively. 

FIG 59 shows a preferred embodiment of a stacked waveguide 
10 phased array applicator having multiple stacked parallel plates 450, 454, 
and 456 forming respective stacked waveguide regions 452, 458, and 
460, each having a corresponding set of monopole phased array antenna 
elements 404 n , 404 n ', 404 n ". 

FIG. 60 shows a block diagram of a preferred embodiment of a 915 
1 5 MHz transmit and receive (T/R) module 800 for use with the monopole 
phased array hyperthermia system 400 (FIG. 47), specifically for use in 
transmit amplifier/phase shift network 1 10 of FIG. 50. Generally, the T/R 
module 800, as well as the monopole hyperthermia system 400, is not 
restricted to operate at 915 MHz, and is adaptable for operation anywhere 
20 within the industrial, scientific, medical (ISM) frequency band of 902 to 
928 MHz. 

The transmit function of the T/R module 800 is used to energize a 
monopole transmit antenna element 404 n (FIG. 50) of the monopole 
hyperthermia array 402 with a 915 MHz signal, having controlled phase 
25 and gain, for the purposes of heating the target 407. Each monopole 
transmit antenna element 404 n of the array 402 is connected to a 
corresponding T/R module and therefore the quantity of T/R modules 
required depends on the quantity of monopole antenna elements 404 in 
the array. n 

30 With regard to the generation of the 91 5 MHz transmit signal, 

having controlled phase and gain, a tunable oscillator 802, tunable from 
180 to 206 MHz, is used to produce a transmit signal having a 193 MHz 
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center frequency. The 193 MHz transmit signal is input to a dual-stage 
voltage-variable attenuator 804 which covers an attenuation range of 0 to 
-40 dB (-60 dB off state) determined by a 12 bit analog control voltage. 
The other port of mixer 808 is driven with a constant frequency 1 108 
5 MHz phase-controlled signal generated by a phase-controlled local 
oscillator 810. 

Phase-controlled local oscillator 810 is excited with a 277 MHz 
signal generated by a fixed frequency local oscillator 812. The output of 
the 277 MHz local oscillator 812 is input to a voltage-variable phase 

10 shifter which shifts the phase of the 277 MHz signal from 0 to 90°. 
determined by a 12 bit analog control voltage. The phase-shifted 277 
MHz signal is input to a X4 frequency multiplier 816 which quadruples the 
signal to 1 108 MHz and extends the phase control range to 0 to 360°. 
The output of the frequency multiplier passes through a 1 108 MHz 

15 bandpass filter 818, having a 50 MHz bandwidth, to remove undesired 
harmonics. The 1 108 MHz output of the bandpass filter is then amplified 
to saturation by an amplifier 820 to produce a relatively constant input 
power to the input port of mixer 808 independent of the commanded 
phase shift. 

20 The upconverted. gain and phase-controlled 91 5 MHz signal output 

from mixer 808 passes through a 915 MHz bandpass filter 822. having a 
50 MHz bandwidth, to remove undesired harmonics generated in mixer 
808. The output signal of bandpass filter 822 passes through another T/R 
switch 824 which, when in the transmit position (T). connects the 
25 bandpass filter output signal to drive the input of a power amplifier 826 
having an average CW output power of up to, or greater than, 100 watts. 
The output signal of the power amplifier 826 passes through another T/R 
switch 828 which, when in the transmit position (T), connects the 
amplified 915 MHz transmit signal to the input of another 915 MHz 
30 bandpass filter 830. having a 50 MHz bandwidth, which removes 

unwanted harmonics generated in the power amplifier. Finally, the 915 
MHz transmit signal output from bandpass filter 830 is connected by a 
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transmission line 832 to a monopole antenna element 404 n of the 
monopole hyperthermia array 402. 

The receive function of the T/R module 800 can be used for passive 
microwave radiometry for non-invasively sensing the temperature of the 
5 target tissue 407 (HQ. 50, with the monopole antenna elements 404 n of 
array 402. To operate in a non-invasive microwave radiometry mode, the 
transmit power is turned off as desired for a period of several seconds 
during which the elements of the monopole hyperthermia array act as 
passive receive antennas. 
10 In this case, the three T/R switches 806, 824, and 828 are set to 

the receive (R, position. The 915 MHz center-frequency passive signal 
received I by the monopole antenna element 404 n is filtered by bandpass 
filter 830, ampiified by a .ow-noise amplifier 834, and again filtered by 

15 ' no*Z 822 " ^ OUtPUt ° f bandP3SS fiter 822 iS mixed the 
15 1 108 MH2 controiled-phase local oscillator signa. by mixer 808, the output 

of wh,ch is input to a lowpass filter 836. Lowpass filter 836 has a 350 
h '' 9h freqUenCy cutoff "«** Prides a 193 MHz center frequency 
receive signa, for input to a second mixer 838. A variabie frequency ,oca, 

20 T^^Tr ' ,50 " 176 ^ ran8e ' a " ° UtPUt -ned 
20 to 163 MHz wh.ch ,s amplified by an amplifier 841 and input to another 

port of mixer 838 to mix with the 193 MHz center frequency receive 

s.gnal. The output of mixer 838 contains a 30 MHz center frequency 

receive signa, which is input to a 30 MHz bandpass filter 842, having a 5 

25 mm to rem ° Ve UOWanted ° Ut ° f band The filtered 30 

25 MHz center frequency receive signal is then passed through a voltage 
controlled 0 to -40 dB attenuator 844 whose output signa, 846 is 
controlled by a 12 bit analog voltage level. 

The 30 MHz center frequency receive signal 846 can be analyzed 
w,th a commercial network analyzer, such as a Hewlett Packard 8510 
analyzer, or can be combined with the output of the other T/R modules 
Ce., receive signals from the other monopole antenna elements, in a 
commercial analog power combiner at the 30 MHz frequency 
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Alternatively, the 30 MHz receive signal 846 can be mixed with a 28.5 
MHz local oscillator and downconverted to a baseband offset frequency of 
1.5 MHz. The resulting baseband signal is lowpass filtered with a cutoff 
frequency of 2.0 MHz, and sampled with a high speed digital to analog 
5 convenor at 4.5 MHz (i.e.. above the Nyquist sampling limit for the 
bandlimited signal). The frequency spectrum of the baseband signal is 
then computed using digital signal processing techniques (ass, j.R. 
Johnson, BLsL. "An Experimental Adaptive Nulling Receiver Utilizing the 
Sample Matrix Inversion Algorithm with Channel Equalization", IEEE 
10 Transactions on Microwave Theory and Techniques, Vol. MTT-39, No. 5, 
pp. 798-808, May 1991). it should further be noted that the T/R module 
800 can be used as a receiver 114 (FIG. 50) for the E-field probes 412 if 
the low-noise amplifier 834 is bypassed. 
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CLAIMS 

A hyperthermia applicator for inducing a temperature rise in a 
target, comprising 

a plurality of electric field radiators; 
a source of electric field radiation coupled to each 
electric field radiator through a controllable transmit weighting 
network coupled to a respective electric field radiator, each 
weighting network controlling the phase and amplitude of the 
electric field radiation transmitted by the respective electric 
field radiator in response to a respective feedback signal; 

at least one electric field probe for detecting electric 
field radiation from the plurality of radiators; and 

a controller coupled to the electric field probes for 
receiving the detected electric field radiation and generating 
the respective feedback signals, and for adjusting the 
feedback signal in response to the detected electric field 
radiation so that the detected electric field radiation is 
minimized at the electric field probe. 

The apparatus of Claim 1, wherein the plurality of electric 
field radiators comprises a phased-array of electric field 
transmit elements. 

The apparatus of Claim 2, wherein the phased-array of 
electric field transmit elements comprises an annular array for 
surrounding the target. 

The apparatus of Claims 1, 2, or 3 wherein the electric field 
probes comprise a plurality of probe elements disposed non- 
invasively along the perimeter of the target. 
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The apparatus of Claims 1, 2, 3, or 4 further comprising 

at least one secondary electric field probe for detecting 
electric field radiation; 

wherein the controller adjusts the feedback signals in 
response to the electric field radiation detected by the 
secondary electric field probe so that the electric field 
radiation is maximized at the secondary electric field probe. 

The apparatus of Claim 5, wherein the plurality of electric 
field radiators comprises a phased-array of electric field 
transmit elements, and the secondary electric field probe is 
disposed at the desired focus of the phased-array. 

The apparatus of Claim 6, wherein 

the secondary electric field probe comprises at least 
one probe element disposed invasively within the target at the 
desired focus of the electric field energy. 

The apparatus of Claims 1, 2, 3, 4, 5, 6 or 7 wherein the 
controller performs a matrix inversion to generate the 
feedback signal. 

The apparatus of Claim 8, wherein the controller comprises 
a receiver coupled to the electric field probes for 

receiving voltages v\,v^ • • .,v£ due to the mth and nth 

electric field radiator at the /th electric field probe; 

means for forming a cross correlation R* of the 

received voltages, the cross correlation being given by 
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where # means complex conjugate and El-) means 
mathematical expectation; 

means for forming a channel correlation matrix R given 

by 

where N aux is the number of electric field probes, /?, is the 
sample correlation matrix observed at the hh electric field 
probe; and 

means for generating an adapted transmit weight 
vector w a given by 

where generally the transmit weight vector 
w=iw 1 ,w 2 ,.„,w N ) T represents the transmit weights of each 
of N controllable weighting networks, and w q represents the 
quiescent transmit weight vector; 

wherein the controller adjusts the feedback signal in 
response to the adapted transmit weight vector w s . 

The apparatus of Claim 1 f 2, 3, 4, 5, 6 or 7 wherein the 
controller performs a gradient search algorithm to generate 
the feedback signal. 



The apparatus of Claim 10, wherein the controller comprises 
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a receiver coupled to the electric field probes for 
receiving voltage amplitudes | v f 'j\, | v 2 J|, . ..,J vfy due to the 
mth and nth electric field radiator at the /th electric field probe 
for the /th configuration of the transmit weights w nJ ; 

means for calculating a figure of merit F nJ from the 
received voltage amplitudes given by 



10 



15 



where N aux is the number of electric field probes, 

means for dithering the transmit weights w nj by a small 

amount in amplitude, AA nJ , and phase, &<t> n/ ; 

means for determining the figure of merit differences 

t^Anj and LF q>nj caused by dithering the amplitude and phase, 
respectively, 

means for determining gradient search directions 
search directions r Anj and r 0 „ y given by 



v- 



AA 



n-1 A A A* 



and 



A* 



«-i A>* A* 
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respectively, and 

means for generating the new transmit weight w n ff+ fJ 
for the (f+ fjth configuration where the amplitude component 
of the new weight is given by 



and the phase component of the new weight is given by 



12. A method for inducing a temperature rise in a target, 
comprising the steps of: 

surrounding a target with an annular phased array of 
electric field radiators; 
10 coupling a source of electric field radiation to each 

electric field radiator through a controllable transmit weighting 
network coupled to a respective electric field radiator; 

controlling the phase and amplitude of the electric field 
radiation transmitted by the respective electric field radiator 
15 with e ach weighting network in response to a respective 

feedback signal- 
detecting electric field radiation from the annular 
phased array of radiators with a plurality of electric field 
probe elements disposed non-invasively along the perimeter 
20 of the target; and 

receiving the detected electric field radiation and 
generating the respective feedback signals for adjusting the 
controllable transmit weighting network in response to the 
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detected electric field radiation so that the detected electric 
field radiation is minimized at the electric field probe 
elements. 



13. A hyperthermia applicator for inducing a temperature rise in a 
target within a body, comprising 

a plurality of electric field radiators; 

a plurality of controllable transmit weighting networks, 
each such network coupled to a respective electric field 
radiator, each weighting network controlling the phase and 
amplitude of the electric field radiation transmitted by the 
respective electric field radiator in response to a respective 
feedback signal; 

a source of electric field radiation coupled to each electric 
field radiator through a respective weighting network; 

at least one electric field probe disposed outside the body 
for detecting electric field radiation from the plurality of 
radiators; and 

a controller coupled to the electric field probes for 
receiving the detected electric field radiation outside the body 
and generating the respective feedback signals, and for 
adjusting the feedback signals in response to the detected 
electric field radiation outside the body so that the electric field 
radiation is controlled at the target inside the body. 

14. The apparatus of Claim 13, wherein the plurality of electric field 
radiators comprises a phased-array of electric field transmit 
elements. 

15. The apparatus of Claim 13 or 14, wherein the electric field 
transmit elements comprises an array of monopole antenna 
elements. 
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16. The apparatus of Claim 15, further comprising an RF reflecting 
groundplane surface for mounting the monopole antenna 
elements, wherein the monopole antenna elements are 
perpendicularly mounted to the same side of the RF reflecting 
ground plane, wherein the body containing the target is 
disposed on the same side of the ground plane as the monopole 
antenna elements. 

17* The apparatus of Claim 16 wherein the ground plane comprises 
an aperture through which the body containing the target is 
disposed to allow positioning the target on the same side of the 
ground plane as the monopole antenna elements. 

18. The apparatus of Claim 16 or 17 further comprising an RF 
reflecting screen mounted perpendicular to the groundplane 
surface behind the monopole antenna elements to reflect RF 
energy from the monopole antenna elements toward the target. 

19. The apparatus of Claim 18 wherein the RF reflecting screen is 
positioned between 1/8 to 1/2 wavelength from the monopole 
antenna elements. 

20. The apparatus of Claim 15, 16, 17, 18 or 19 further comprising 
an enclosure surrounding the monopole antenna elements and 
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providing a vessel for enclosing a bolus of fluid between the 
monopole antenna elements and the body. 

21 . The apparatus of Claim 20 wherein the vessel comprises a 
bolus of deionized water. 

22. The apparatus of Claim 15, 16, 17, 18, 19, 20, or 21 wherein 
the monopole antenna elements are arranged in a circular arc of 
substantially constant radius. 

23. The apparatus of Claim 22 wherein the body containing the 
target is a cranium and the radius is substantially the distance 
from the monopole antenna array to the center of the cranium. 

24. The apparatus of Claim 22 wherein the body containing the 
target is a cranium and the radius is substantially the distance 
from the monopole antenna array to the target within the 
cranium. 

25. The apparatus of Claim 22, 23, or 24 wherein the monopole 
antenna elements resonate at between 800 and 1000 MHz. 

26. The apparatus of Claim 13, 14, 15, 16, 17, 18, 19, 20, 21, 
22, 23, 24, or 25 wherein the electric field probes comprise a 
plurality of probe elements disposed non-invasively along the 
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perimeter of the body between the elements of the phased- 
array and the target. 

27. The apparatus of Claim 26 wherein the controller comprises 
means for adjusting the feedback signals to minimize the 
difference in the electric field detected by adjacent electric field 
probes and thereby provide uniform electric field radiation into 
the body. 

28. The apparatus of Claim 27, wherein the adjusting means 
performs a gradient search algorithm to adjust the feedback signals. 

29. The apparatus of Claim 27 f wherein the adjusting means 
performs a matrix inversion algorithm to adjust the feedback 
signals. 

30. The apparatus of Claim 26, wherein 

the controller comprises means for adjusting the 
feedback signals to provide a particular electric field pattern 
across the electric field probe array and thereby focus radiation 
into the target. 
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